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Asstract.—Asplenium serratifolium (Aspleniaceae), a new fern species from central Vietnam, is 
described and illustrated. The new species is characterized by plants 10-18 cm tall, laminae 
pinnatipartite, lobe margins entire or with shallow teeth, and veins simple or forked. Molecular 
phylogenetic analysis based on five plastid markers (atpB, rbcL, rps4, rps4-trnS, and trnL-F) 
indicate that the new species is closely related to A. ensiforme. 


Key Worps.—Asplenium ensiforme, pinnatipartite, fern phylogeny, Quang Tri Province 


The largest fern genus, Asplenium L. (Aspleniaceae), is sub-cosmopolitan 
and comprises more than 700 species (Lin and Viane, 2013), of which 44 
species are known to occur in Vietnam (Pham-Hoang, 1999). During fieldwork 
to survey ferns and lycophytes in central Vietnam in November 2014, two 
collections of a species of Asplenium were obtained from Quang Tri Province. 
The pinnatipartite lamina of this species distinctly differs from that of any 
other Asplenium known in Asia (Pham-Hoang, 1999; Lin and Viane, 2013; 
Lindsay and Middleton, 2012). Based on the investigation of morphological 
characters, literature, specimen comparisons, and phylogenetic analysis, we 
found that this species of Asplenium is new to science (Pham-Hoang, 1999; Lin 
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and Viane, 2013; Lindsay and Middleton, 2012). We here describe this new 
species as Asplenium serratifolium. 


MATERIAL AND METHODS 


Materials and taxon sampling.—We sequenced one individual from each of 
the two known populations of the new species, both from the Bac Huong Hoa 
Nature Reserve, Quang Tri, Vietnam. A total of 26 additional species of 
Asplenium were included in the molecular analyses based on morphological 
similarities and BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PAGE_ 
TYPE=BlastSearch) search results of the rbcL sequence of the new species 
with default settings in BLAST. Three important references, Schneider (2004), 
Schneider et al. (2017), and Ohlsen et al. (2015), were also consulted in 
choosing the ingroup taxa. Species of Hymenasplenium Hayata were used as 
outgroups because there are only two genera in the family generally recognized 
(Ohlsen et al., 2015). 

Morphological study.—The morphological description of the new species 
was based on the holotype and a paratype deposited in the Missouri Botanical 
Garden herbarium (MO). Measurement of the characteristics of roots, stipes, 
rachises, scales, and indusia was conducted using a micrometer under a 
dissecting microscope. Herbarium investigations were carried out at CDBI, 
MO, PHH, SYS, and VNMN and by using online images in the Chinese Virtual 
Herbarium (http://www.cvh.org.cn) provided by various Chinese herbaria and 
on those in JSTOR Global Plants (https://plants.jstor.org/) provided by herbaria 
worldwide. | 

DNA extraction, PCR amplification and sequencing.—Total genomic DNA 
was extracted from silica-gel-dried leaves using the modified 2XCTAB 
procedure of Doyle and Doyle (1987). DNA sequence data were obtained 
from five plastid markers: atpB, rbcL, rps4, rps4-trnS intergenic spacer, and 
trnL-F intergenic spacer. Primer sequences and PCR procedures followed Xu et 
al, (2018). Amplified fragments were purified with TIANquick Mini 
Purification Kits (TIANGEN). Purified PCR products were sequenced by 
TSINGKE Biological Technology (Guangzhou, China). DNA sequences of 23 
Asplenium species were downloaded from GenBank. All sequences used in 
this study, together with their GenBank accession numbers, are listed in 
Appendix I. 

Sequence alignment and phylogenetic analysis.—Sequences were edited 
and assembled using Sequencher V.4.14 (GeneCodes Corporation, Ann Arbor, 
Michigan, USA). Sequence alignment and adjustment were conducted 
manually with Bioedit (Hall, 1999). Phylogenetic relationships were 
analyzed through maximum parsimony (MP), maximum likelihood (ML), 
and Bayesian inference (BI). MP tree searches were carried out for each data 
matrix in PAUP v. 4.0b10 (Swofford, 2002). Parsimony jackknife (JK) analyses 
(Farris et al., 1996) were conducted using PAUP with the removal probability 
set to approximately 37%, and “jac” resampling emulated. Two hundred 
replicates were performed from different random addition sequence trees, with 
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TaBLE 1. Best-fitting models, data matrices, and tree statistics for each of the analyses. Missing 
data include missing sequences, uncertain bases (N, R, Y, V, etc.), and gaps (-). 


Locus | Selected model # accessions # missing (%)  # chars. 
Plastid atpB gene TIM1+G 6 84.5 1239 
Plastid rbcL gene TIM1+G 37 2.1 1252 
Plastid rps4 gene & rps4-trnS spacer TPM1uf+I+G 18 63.5 1030 
Plastid trnL intron & trnL-F spacer TVM+G 18 62.4 384 
Simultaneous TIM1+I+G 37 50.4 3905 


100 TBR searches per replicate and a maximum of 100 trees held per TBR 
search. ML tree searches and ML bootstrapping (BS) were performed using the 
web server RAXML-HPC2 on XSEDE v.8.2.10 on the CIPRES web server with 
1,000 bootstrap replicates (Miller, Pfeiffer, and Schwartz, 2010). Bayesian 
inference was conducted using MrBayes 3.1.2 (Ronquist and Huelsenbeck, 
2003) on CIPRES (Miller, Pfeiffer, and Schwartz, 2010). Two runs of four 
Markov chain Monte Carlo chains were conducted, each beginning with a 
random tree and sampling one tree every 1000 generations of 10,000,000 
generations. Convergence among chains was checked using Tracer ver. 1.4 
(Rambaut and Drummond, 2007), and the first 25% was discarded as burn-in to 
ensure that stationarity had been reached. The remaining trees were used to 
calculate a 50% majority-rule consensus topology and posterior probabilities 
(PP). 

jModelTest 2 (Guindon and Gascuel, 2003; Posada, 2008; Darriba et al., 
2012) was used to select the best-fitting likelihood model. The Akaike 
information criterion (Akaike, 1974) was used to select among models (Table 
1), following Pol (2004) and Posada and Buckley (2004). 


RESULTS AND DISCUSSION 


The alignment of the combined plastid DNA datasets (atpB: 1,239 bp, rbcL: 
1,252 bp, rps4 and rps4-trnS intergenic spacer: 1,030 bp, trnL-F intergenic 
spacer: 384 bp) was 3,905 characters long (details concerning the datasets 
analyzed and statistics for the resulting trees are shown in Table 1) of which 
535 were parsimony-informative and 388 were variable and parsimony- 
uninformative. In our plastid phylogeny, the three phylogenetic analyses (BI, 
ML, and MP) revealed similar topologies, with minor differences at some 
weakly supported branches, and all three trees indicated that the new species 
formed a strongly supported clade (MLBS: 100%, BIPP: 1, MPJK: 100%). 

The two samples of the new species had identical DNA sequences, which 
were different from all other sampled species. In addition, the phylogenetic 
analyses indicate that the new species is closely related to A. ensiforme 
Wallich ex Hooker & Greville (Hooker and Greville, 1828; Fig. 1). Morpholog- 
ically, A. serratifolium can be easily distinguished from A. ensiforme in having 
plants 10-18 cm tall, laminae pinnatipartite, and lobe margins entire or with 
shallow teeth, whereas A. ensiforme has plants 30-45 (—65) cm tall, laminae 
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Asplenium tenuicaule TNS: 765223 Yamagata, Japan 
Asplenium scortechinii Fan 01144 Hainan, China 

90/1/100.| Asp/enium sarelii TNS: 769194 Mie, Japan 
400/1/100 | 'Asplenium sarelii Nogami & Oohora 11 Mie, Japan 
Asplenium pekinense Viane 10031 Beijing, China 
Asplenium pekinense Lu SG/C67 Yunnan, China 
Asplenium varians Fraser-Jenkins 10046-10047 China 
Asplenium simaoense XKW317-1 Yunnan, China Clade X 
Asplenium gueinzianum Chengs.n. China 
Asplenium tenuifolium JLO0361 Yunnan, China 
400/1/100 | ASPlenium yunnanense_ Fraser-Jenkins 10044-10045 China 
Asplenium yunnanense Lu SG/D22 Yunnan, China 
Asplenium lushanensé Lu SG/D21 Yunan, China 
Asplenium interjectum Lu SG/J12 Yunnan, China 
100/1/100 ae aires pulcherrimum TTJ-GZ-015 Yunan, China 4 
Asplenium coenobiale TNS: 763925 Kochi, Japan 
99/1/100 | AsPlenium ensiforme Fukuoka et al. 93-T619 Doi Inthanon, Thailand 
Asplenium ensiforme Nogami etal. 8 Mie, Japan 
Asplenium ensiforme TNS: 763995 Miyazaki, Japan 
Asplenium serratifolium sp nov Zhang et al. 7631 Quang Tri, Vietnam 
100/1/100 ' Asplenium serratifolium sp nov Zhang et al. 7614 Quang Tri, Vietnam 
Asplenium yoshinagae Lu SG/A4 Yunnan, China 
Asplenium pseudowilfordii Nogami 14 Nara, Japan 

100/1/100 + Asplenium wilfordii TNS: 763478 Kagoshima, Japan 

Asplenium wrightii Cranfill TW040 Taiwan Island 
Asplenium wrightioides Lu SG/J13 Yunnan, China | Clade VII 
Asplenium ceterach Van De Riet #707 Unknown | Clade IV 
Asplenium antiquum P22088 Unknown 
Asplenium nidus Berkeley 680392 Madagascar 
Asplenium griffithianum Fujita Yaku-1204 Japan 

| Asplenium griffithianum Murakami J93-001 Japan 
99/1/100 | Asplenium griffithianum Murakami J93-001 Japan 
Asplenium prolongatum Cult. in NYBG 
Hymenasplenium unilaterale Hemp A. 18 Kenya 
400/1/100 Hymenasplenium excisum Brownsey & Perrie FiJi 190 Fiji 

86/-/100 Hymenasplenium cheilosorum Schater 55 Yunnan, China 
97/1/1400 Hymenasplenium obliquissimum Murakami & Cheng 94-M 1811 Yunnan, China 


100/1/100 


95/1/80 


100/1/100 


100/1/100 


100/1/100 


100/1/100 


New Clade 
100/1/100 


100/1/100 


-/-194 
100/1/100 


-/0.94/- 


100/1/100 


87/0.99/100 | Thamnopteris 1 


100/1/100 


100/1/100 


A. tenerum clade 
400/1/100 


Outgroups 
“Ne 


0.02 


Fic. 1. The phylogenetic position of Asplenium serratifolium based on five plastid markers (atpB, 
rbcL, rps4, rps4-trnS, trnL-F). The numbers associated with branches are maximum likelihood 
bootstrap support (MLBS), Bayesian posterior probability (BIPP), followed by maximum parsimony 
jackknife support (MPJK). “-” indicates MLBS < 85%, BIPP < 0.9%, or MPJK < 80%. The major 
clades (either new or following Ohlsen (2015)) are indicated on the right. 


simple and entire to repand. Detailed comparison of the two species is listed in 
Table 2. 


Only two populations of the new species in the Bac Huong Hoa Nature 


Reserve of Vietnam were found, while the related species Asplenium 
ensiforme is widely distributed in Asia (Lin and Viane, 2013). 


TaBLeE 2. Morphological comparison between Asplenium serratifolium and A. ensiforme. 


Character A. serratifolium A. ensiforme 
Plant height 10-18 cm 30—45(—65) cm 
Petiole length 1-6 cm 5—-8(—15) cm 
Rhizome scales 2-4 X 0.08-0.15 mm Ca. 6 X 1mm 
Lamina Pinnatipartite - Simple 
Lobe margins Entire or serrated Entire to repand 
Distribution Central Vietnam Southeast Asia 


Veins Simple or forked Once-forked 
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TAXONOMIC TREATMENT 


Asplenium serratifolium Li Bing Zhang & K.W. Xu, sp. nov. (Figs. 2, 3). TYPE: 
VIETNAM. Quang Tri: Huong Hoa District, Bac Huong Hoa Nature 
Reserve, Huong Viet Commune, Ta Rung Village, 1250 m, 16°48’N, 
106°35’E, 22 Nov 2014, Li-Bing Zhang, Liang Zhang & Ngan Thi Lu 7631 
(holotype: MO!; isotypes: CDBI!, VNMN!}). 


Diagnosis.—Asplenium serratifolium somewhat resembles A. ensiforme by 
its stipes sulcate adaxially, midribs raised abaxially and sulcate adaxially, 
obscure veins, and sori in rows on two sides of midrib separately. The former 
has plants 10-18 cm tall, laminae pinnatipartite, lobe margins entire or with 
shallow teeth, and veins simple or forked, while the latter has plants 30—45(- 
65) cm tall, laminae simple, entire to repand, and veins once-forked. 

Plants perennial, evergreen, 10-18 cm tall. Rhizomes erect, short, apices 
densely scaly; scales brown to dark brown, lanceolate, 2-6 X 0.18-0.5 mm, 
margins entire. Fronds simple and caespitose, herbaceous, yellowish green or 
grayish green when dry, subglabrous; stipes stramineous, sulcate adaxially, 
rounded to carinate abaxially, 1-6 cm, bases densely covered with scales, 
becoming subglabrous above, scales clathrate, brown, lanceolate, 2-4 x 0.08— 
0.15 mm, margins entire; laminae pinnatipartite, lanceolate, 5-12 < 0.6—1.8 
cm, bases reduced and decurrent on stipes, apices acute to acuminate, widest 
at middle part of laminae, margins of lobes entire or serrate, teeth shallow and 
obtuse. Midribs stramineous, adaxially sulcate, raised abaxially, sparsely 
scaly, scales brown, lanceolate, ca. 0.08 X 0.8 mm, entire, glabrous; veins 
obscure, simple or forked. Sori linear, (0.1—) 0.2-0.5 cm, starting close to 
midribs at an angle of 25°—45° on acroscopic veinlets; indusia yellowish green 
or stramineous and becoming brown when dry, linear, papery, margins entire, 
opening toward midribs, persistent. 

Geographical distribution.—Asplenium serratifolium is known from two 
locations within the Bac Huong Hoa Nature Reserve, Quang Tri, central 
Vietnam. 

IUCN Red List category.—Only two populations and fewer than 50 
individuals of Asplenium serratifolium are known. We did not find any 
other populations in the course of our three expeditions (from 2013-2015) to 
collect Vietnamese pteridophytes. As a result, a CR-Critically Endangered 
category is suggested based on current information available and following 
IUCN (the International Union for Conservation of Nature and Natural 
Resources) guidelines (IUCN, 2015). 

Etymology.—Based on the Latin prefix, serrati-, serrate, and the Latin suffix, 
-folium, leaf, referring to the saw-toothed laminae of the new species. 

Additional specimens examined (paratypes).—VIETNAM. Quang Tri: 
Huong Hoa District, Bac Huong Hoa Nature Reserve, Huong Viet Commune, 
Ta Rung Villige, Sa Mu Pass, elev. 1300 m, 16°48’N, 106°35’E, 21 Nov 2014, Li- 
Bing Zhang, Liang Zhang & Ngan Thi Lu 7614 (CDBI!, MO!, VNMN!). 
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Fic, 2. Asplenium serratifolium sp. nov. A and B. Habit. C. Plant. D. Abaxial lamina. E. Sulcate 
stipe, adaxial view. F. Stipe scales. G. Adaxial lamina showing the sulcate midrib. H. Portion of 
abaxial lamina; red arrow shows the obscure, forked veins. 
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Fic. 3. Asplenium serratifolium sp. nov. A. Habit. B. Lamina, abaxial surface, showing simple and 
forked veins. (drawings by Lei Jiang, based on holotype at MO). 
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APPENDIX I. SPECIES USED IN MOLECULAR PHYLOGENETIC ANALYSES, THEIR VOUCHERS, 
AND GENBANK ACCESSION NUMBERS. “—’ MEANS MISSING DATA. 


Asplenium antiquum Makino, P22088 (WELT), Unknown: rbcL EU240033, 
atpB —, rps4 & rps4-trnS EU240020, trnL-F EU240028 (Shepherd et al., 2008). 
Asplenium ceterach L., Van De Riet 707 (EB), Unknown: rbcL HQ676494, 
atpB —, rps4 & rps4-trnS —, trnL-F HQ676516 (de Groot et al., 2011). 
Asplenium coenobiale Hance, TNS: 763925 (TNS), Kochi, Japan: rbcL 
AB574855, atpB —, rps4 & rps4-trnS —, trnL-F — (Ebihara et al., 2010). 
Asplenium ensiforme Wall. ex Hook. & Grev., TNS: 763995 (TNS), Miyazaki, 
Japan: rbcL AB574856, atpB —, rps4 & rps4-trnS —, trnL-F — (Ebihara et al., 
2010). Asplenium ensiforme Wall. ex Hook. & Grev., Nogami et al. 8 (KYO), 
Mie, Japan: rbcL AB014689, atpB —, rps4 & rps4-trnS —, trnL-F — (Murakami 
et al., 1999). Asplenium ensiforme Wall. ex Hook. & Grev., Fukuoka et al. 93- 
T619 (KYO), Doi Inthanon, Thailand: rbcL AB014709, atpB —, rps4 & rps4- 
trnS —, trnL-F — (Murakami et al., 1999). Asplenium griffithianum Hook., 
Fujita et al., Yaku-1204 (FU), Kagoshima, Japan: rbcL AB574857, atpB —, rps4 
& rps4-trnS —, trnL-F — (Ebihara et al., 2010). Asplenium griffithianum 
Hook., Murakami J93-001 (TI), Kagoshima, Japan: rbcL AB013252, atpB —, 
rps4 & rps4-trnS —, trnL-F — (Murakami et al., 1999a). Asplenium 
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griffithianum Hook., Murakami J93-001 (KYO), Kagoshima, Japan: rbcL 
AB014688, atpB —, rps4 & rps4-trnS —, trnL-F — (Murakami et al., 1999a). 
Asplenium interjectum Christ, Lu SG/J 12 (PYU), Yunnan, Ching: uce 
AY545480, atpB —, rps4 & rps4-trnS AY725043, trnL-F AY725038 (Li and Lu, 
2006). Asplenium IJaciniatum D. Don, Cheng s.n. (BM), China: rbcL AY549747, 
atpB —, rps4 & rps4-trnS AY549801, trnL-F AY549851 (Schneider et al., 2005). 
Asplenium lushanense C. Chr., Lu SG/D 21 (PYU), Yunnan, China: rbcL 
AY545481, atpB —, rps4 & rps4-trnS AY725042, trnL-F AY725033 (Li and Lu, 
2006). Asplenium nidus L., U.C. Botanical Garden 68.0392 (UC), Madagascar: 
rbcL AF525270, atpB —, rps4 & rps4-trnS AY549807, trnL-F AF525246-(Pinter 
et al., 2002). Asplenium pekinense Hance, Viane 10031, Unknown: rbcL 
GU929864, atpB —, rps4 & rps4-trnS —, trnL-F — (Unpublished data). 
Asplenium pekinense Hance, Lu SG/C 67 (PYU), Yunnan, China: rbcL 
AY545479, atpB —, rps4 & rps4-trnS AY725040, trnL-F AY725037 (Li and 
Lu, 2006). Asplenium prolongatum Hook., Cult. in NYBG (UC), New York, US: 
rbcL AY549752, atpB —, rps4 & rps4-trnS AY549813, trnL-F AY549856 
(Schneider et al., 2005). Asplenium pseudowilfordii Tagawa, Nogami 14 
(KYO), Nara, Japan: rbcL AB014696, atpB —, rps4 & rps4-trnS —, trnL-F — 
(Murakami et al., 1999). Asplenium pulcherrimum (Baker) Ching ex Tardieu, 
Ke-Wang Xu TTJ-GZ-015 (SYS), Guizhou, China: rbcL MH109142, atpB 
MH109136, rps4 & rps4-trnS MH109148, trnL-F (Unpublished data). Asple- 
nium sarelii Hook., TNS:769194 (TNS), Mie, Japan: rbcL AB574873, atpB —, 
rps4 & rps4-trnS —, trnL-F — (Ebihara et al., 2010). Asplenium sarelii Hook., 
Nogami & Oohora 11, Mie, Japan: rbcL AB014693, atpB —, rps4 & rps4-trnS —, 
trnL-F — (Murakami et al., 1999). Asplenium scortechinii Bedd.; Fan 01144, 
Hainan, China: rbcL MH509434, atpB —, rps4 & rps4-trnS —, trnL-F — (this 
study). Asplenium serratifolium Li Bing Zhang & Ke-Wang Xu, Zhang et al. 
7631 (CDBI, MO, VNMN), Quang Tri, Vietnam: rbcL MH509432, atpB 
MH509430, rps4 & rps4-trnS MH509435, trnL-F MH509437 (this study). 
Asplenium serratifolium Li Bing Zhang & Ke-Wang Xu, Zhang et al. 7614 
(CDBI, MO, VNMN), Quang Tri, Vietnam: rbcL MH509433, atpB MH509431, 
rps4 & rps4-trnS MH509436, trnL-F MH509438 (this study). Asplenium 
simaoense Ke-Wang Xu, Li Bing Zhang & W. B. Liao, Ke-Wang Xu XKW317- 
1 (SYS), Yunnan, China: rbcL MH109137, atpB MH109131, rps4 & rps4-trnS 
MH109143, trnL-F MH109149 (Unpublished data). Asplenium tenuicaule 
Hayata, TNS: 765223 (TNS), Yamagata, Japan: rbcL AB574878, atpB —, rps4 & 
rps4-trnS —, trnL-F — (Ebihara et al., 2010). Asplenium tenuifolium D. Don, 
Lei Jiang JLO0361 (SYS), Yunnan, China: rbcL MH109141, atpB MH109135, 
rps4 & rps4-trnS MH109147, trnL-F MH109153 (Unpublished data). Asple- 
nium varians Wall. ex Hook. & Grev., Fraser-Jenkins 10046-10047 (BM), 
China: rbcL AY300147, atpB —, rps4 & rps4-trnS AY549802, trnL-F AY300094 
(Schneider et al., 2004). Asplenium wilfordii Mett. ex Kuhn, TNS: 763478 
(TNS), Kagoshima, Japan: rbcL AB574883, atpB —, rps4 & rps4-trnS —, trnL-F 
— (Ebihara et al., 2010). Asplenium wrightii D. C. Eaton ex Hook., Cranfill 
TW040 (UC), Taiwan: rbcL AY549730, atpB —, rps4 & rps4-trnS AY549766, 
trnL-F AY549833 (Schneider et al., 2005). Asplenium wrightioides Christ, Lu 
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SG/J13 (PYU), Yunnan, China: rbcL AY725031, atpB —, rps4 & rps4-trnS 
AY725044, trnL-F AY725032 (Li and Lu, 2006). Asplenium yoshinagae 
Makino, Lu SG/A4 (PYU), Yunnan, China: rbcL AY725030, atpB —, rps4 & 
rps4-trnS AY725045, trnL-F AY725036 (Li and Lu, 2006). Asplenium 
yunnanense Franch., Lu SG/D 22 (PYU), Yunnan, China: rbcL AY545482, 
atpB —, rps4 & rps4-trnS AY725041, trnL-F AY725034 (Li and Lu, 2006). 
Asplenium yunnanense Franch., Fraser-Jenkins 10044-10045 (BM), China: 
rbcL AY300149, atpB —, rps4 & rps4-trnS AY549803, trnL-F AY300096 
(Schneider et al., 2004). Hymenasplenium cheilosorum (Kunze ex Mett.) 
Tagawa, Schater 55 (GOET), Yunnan, China: rbcL JF832071, atpB —, rps4 & 
rps4-trnS —, trnL-F — (Rothfels et al., 2012). Hymenasplenium excisum (C. 
Presl) S. Lindsay, Brownsey & Perrie FIJI 190 (WELT), Fiji: rbcL KP774884, 
atpB —, rps4 & rps4-trnS —, trnL-F — (Ohlsen et al., 2015). Hymenasplenium 
obliquissimum (Hayata) Sugimoto, Murakami & Cheng 94-M 1811, Yunnan, 
China: rbcL AB016186, atpB —, rps4 & rps4-trnS —, trnL-F — (Murakami et al., 
1998). Hymenasplenium unilaterale (Lam.) Hayata, Hemp 18 (BM), Kenya: 
rbcL AF240652, atpB —, rps4 & rps4-trnS —, trnL-F — (Pinter et al., 2002). 
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AsstTRACT.—The transition zone from the tropics to the subtropics in Mexico is an outstanding and 
poorly explored area of biotic overlap, and as such very interesting to evaluate how species richness 
patterns change in the latitudinal context. In order to assess alpha diversity patterns of ferns in this 
zone, we analyzed changes in their assemblages along an elevational gradient from 300 m to 3,000 m 
in the Sierra de Juarez, Oaxaca, and compared the results with similar studies conducted in Central 
and South America. We report the most detailed plot-based elevational gradient study on ferns in the 
region and provide detailed data on species abundances and elevational ranges. In 121 plots of 20 x 
20 m, we found 195 fern species, with the highest species numbers recorded at mid-elevations 
(1,100—2,200 m) and lower richness at both gradient extremes. This hump-shaped pattern was 
mainly driven by epiphytes, which contributed about 40% to the total species richness and had 
wider elevational ranges than terrestrial species, which contributed 38% to the total richness, and 
were more or less constant at low to mid elevations about 2,500 m, followed by a gradual decline. 
Overall, richness per plot was low compared to other elevational gradients within the tropics closer 
to the Equator. This decline of species richness towards the subtropics and temperate regions 
appears to be rather abrupt within Central America, but additional data are needed to close the 
remaining gap of knowledge between Costa Rica and Mexico. 


Key Worps.—Pteridophytes, alpha diversity, distribution, elevational range, Priority Terrestrial 
Region. 
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Ferns and lycophytes (although different lineages, henceforth called ferns) 
globally comprise about 11,000 species, thus being the second most diverse 
group of vascular plants on Earth (Moran, 1998; PPG I, 2016). With more than 
1,000 recorded species, they are well represented in Mexico, which has one of 
the best-documented fern floras in the world (Mickel and Smith, 2004). The 
country is considered an epicenter for species radiations in some fern groups 
(principally groups specialized in seasonally dry habitats), and the states of 
Oaxaca, Chiapas, and Veracruz harbor the highest species numbers in the 
country (Krémer et al., 2015; Mickel and Beitel, 1988; Mickel and Smith, 
2004). Because fern taxonomy is well studied and non-woody life forms 
contribute up to 58% of the vascular plant richness in tropical forests (Gentry 
and Dodson, 1987; Linares-Palomino et al., 2009), ferns represent a valuable 
group for ecological and biogeographic studies. 

In this regard, elevational gradients are an interesting natural experiment, 
because they include major environmental changes within short horizontal 
distances, while dispersal limitation is minimized. Ferns have been widely 
studied along tropical elevational gradients in recent years, and show 
consistent general richness pattern at local and regional scales, with highest 
richness commonly found at mid-elevations and decreasing diversity towards 
both high and low elevations (Bhattarai, Vetaas, and Grytnes, 2004; Karger et 
al., 2011; Kessler, 2001; Kessler et al., 2011; Kluge, Kessler, and Dunn, 2006; 
Parris, Beaman, and Beaman, 1992; Salazar et al., 2013; Watkins et al., 2006). 

With respect to major life-forms, this hump-shaped distributional trend in 
tropical mountains is commonly driven by epiphytes, which are extraordi- 
narily rich in multi-structured tree crowns within the humid cloud 
condensation layer between 1,000—3,000 m (Acebey, Krémer, and Kessler, 
2017; Kessler, 2000; Kluge and Kessler, 2005; Kromer et al., 2013; Watkins et 
al., 2006). 

However, how elevational fern richness changes with latitude is less well 
explored, and two alternative hypotheses have been proposed. The first 
assumes that richness peaks along elevational gradients gradually shift 
towards lower elevations with increasing distance from the Equator in parallel 
with decreasing temperatures. Thus, Kessler et al. (2011) showed that the 
maximum richness of ferns along elevational gradients in tropical wet 
environments is situated at around 15-17°C mean annual temperature, and 
this level should be found at lower elevations at higher latitudes. Alternative- 
ly, when comparing various fern transects from Bolivia to Mexico, Salazar et 
al. (2013) found that overall elevational richness is strikingly constant between 
Costa Rica and northern Bolivia. The major shift towards lower species 
richness was found at the transition from tropical to subtropical climates. 
However, this decline was not linked to a downward shift in the elevation of 
the richness peak, but rather reflected a general decline of species richness, as 
also found by Krémer et al. (2013), Acebey, Kromer, and Kessler (2017), and 
Carvajal-Hernandez et al. (2017) in the regions of Los Tuxtlas and central 
Veracruz, in the southeast of Mexico. Apart from these publications, no other 
suitable transects have been studied thus far to allow a fine-scale assessment of 
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Fic. 1. Study area in the Sierra de Juarez, Oaxaca, Mexico, showing the locations of the sampling 
sites (dots). 


the change of species richness patterns at the Mexican transition zone (MTZ, 
sensu Halffter and Morrone, 2017) between tropical and subtropical climates. 

In the present study, we attempt to close a gap in our knowledge about 
elevational species richness patterns of ferns starting with the analysis of an 
extensive elevational gradient in the Sierra de Juarez, Oaxaca, inside the MTZ 
by presenting plot-based species richness data together with an annotated 
species list with detailed information on abundance and life form assignments. 
Specifically, we ask the following questions: 


1. How does local species richness of ferns change with elevation? 

2. Is the elevational pattern of species richness and ranges different between 
both major life forms (terrestrials and epiphytes)? 

3. Is the total and maximum species richness significantly lower in magnitude 
and elevation compared to other studies performed in Central and South © 
America? 


METHODS 


Study area.—The study area was located in southeastern Mexico in the 
Sierra de Juarez (also called Sierra Norte), Oaxaca, and the majority of the 
localities belong to the Chinantla region (Fig. 1, additional plot information in 
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Appendix 1). The Sierra de Juarez has been declared a Priority Terrestrial 
Region due to its species and ecosystem richness (Arriaga et al., 2000), and 
thus is an important area for conservation. The eastern slope of the Sierra de 
Juarez hosts a very rich and one of the most distinctive floras in Mexico (Meave 
et al., 2017) and represents the northernmost location within the Neotropics, 
where a large and well-preserved elevational gradient of humid tropical forest 
still remains (Boyle, 1996). There is an abrupt elevational increase in this area 
in less than 30 km from Valle Nacional (ca. 65 m) to the top of Humo Chico 
Mountain (ca. 3,200 m; Meave, Rincén, and Romero-Romero, 2006), called 
the—fern hunter‘s paradise by John Mickel (1965) during fieldwork to prepare 
the first treatment about ferns from Oaxaca (Mickel and Beitel 1988). Given the 
extensive deforestation elsewhere in Mexico (Rosete-Vergés et al., 2014), the 
presence of this gradient of humid montane forest in Oaxaca is remarkable 
(Rzedowski and Palacios-Chavez, 1977). Despite a generally good conservation 
status above 1,000 m, human activities are leaving an important footprint in 
the lower parts of the area (Meave et al., 2017). In addition, the Sierra is 
periodically impacted by strong winds from the Gulf coast (>80 km/h), which 
disturb the forest structure (Alfonso-Corrado et al., 2017). There are few 
general studies on vegetation along the whole gradient in the Sierra de Juarez 
(Boyle, 1996; Meave, Rincén, and Romero-Romero, 2006; Rzedowski and 
Chavez, 1977; Rincdén-Gutierrez, 2007). However, a recently published 
extensive floristic survey in the region reported 130 fern species along an 
elevational range from 250—3,020 m (Meave et al., 2017). Forest zonation (with 
variation) begins in the lowlands with tropical rain forest from 300—1,000 m, 
followed by a transitional forest zone from 1,100—1,300 m, and cloud forest at 
1,400—2,000 m (sometimes to 2,500 or more). Within the cloud forest belt, there 
are areas of forest dominated by Oreomunnea mexicana. Above 2,500 m there 
is Quercus forest and above 3,000 m elfin forest forms the treeline ecotone at 
around 3,100 m. 

Climate.—La Chinantla is considered a hyperhumid region in Mexico 
(Meave, Rincon, and Romero-Romero, 2006, Meave et al., 2017), but there are 
few climatic records for localities in the Sierra de Juarez. The climate at the 
lowest elevation is hot (>22°C) changing to a warm climate (18—22°C) higher 
up, and temperate at higher elevations (5-12 °C) (Meave, Rincon, and Romero- 
Romero, 2006). Mean annual precipitation is around 3,590 mm at the lower 
part of the gradient (Valle Nacional), 5,800 mm in the middle (La Esperanza), 
and 3,240 mm at the upper part at the locality of Humo Chico (Meave, Rincén, 
and Romero-Romero, 2006). Due to the scarcity of meteorological records in 
the study area, we extracted ecologically important temperature and 
precipitation data (Fig. 2) from the global climate dataset CHELSA (Karger et 
al., 2017). 

Fern sampling.—A total of 121 plots were inventoried at 15 sites along an 
elevational gradient between 300 m and 3,000 m at elevational intervals of 
200-400 m (Appendix 1). At each elevation, eight plots of 20 X 20 m were 
sampled with a consistent, standardized methodology (Karger et al., 2014; 
Kessler and Bach, 1999). The plots were established in natural zonal forest, 
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Fic. 2. Elevational trends of (a) mean annual temperature and (b) precipitation. Data derived from 
the CHELSA data set (Karger et al., 2017) for each plot of the elevational transect in the Sierra de 
Juarez, Oaxaca, Mexico. Trend lines were set with locally weighted regressions (LOWESS). 


avoiding special structural features like canopy gaps, ridges, ravines, riparian 
areas, treefall gaps, landslides, and other disturbed areas whenever possible. In 
each plot, all fern species and their abundances were recorded for terrestrials 
(growing on soil, rocks, and dead wood) and epiphytes. In most cases, the 
abundance (i.e., number of individuals of a species) is easy to assess, 
exceptions are only small species with distant fronds on long creeping and 
interwoven rhizomes (e.g., Hymenophyllaceae). Because in these cases exact 
numbers of individuals can only be obtained by destructive sampling, we 
treated a species patch as one individual in order not to overestimate the 
abundance of these species. Epiphytes were sampled up to heights of eight m, 
with trimming poles and recorded using binoculars at greater heights, by 
climbing lower parts of trees, and by searching recently fallen trees and 
branches within and adjacent to the plots (Gradstein et al., 2003). 

Samples of all fern species were collected and deposited in the University 
Herbarium, University of California (UC) in Berkeley, USA, herbarium XAL of 
the Instituto de Ecologia, A. C. (Xalapa, Veracruz, Mexico) and herbarium 
MEXU of the Universidad Aut6noma de México (Mexico City, Mexico). 
Collections were identified by the first author and Alan R. Smith (UC). 
Taxonomy primarily followed Mickel and Smith (2004) and the current 
classification for lycophytes and ferns established by the Pteridophyte 
Phylogeny Group (PPG I, 2016) and Hassler (2018). 

Data analysis.—We assessed plot-based alpha-diversity for all species and 
separately for the two major lifeforms (terrestrials and epiphytes) and 
described their trends by locally weighted regressions (LOWESS). In order to 
assess whether the elevational ranges of terrestrials and epiphytes were 
different, we used a Wilcoxon signed-rank test. We considered a species as 
belonging to both life forms when more than 20% of the total number of 
individuals occurred in each life form. The frequencies of these three 
categories of life forms were compared with an analysis of variance (ANOVA) 
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Fic. 3. Richness patterns of ferns (number of species per plot) along the elevational gradient in 
Sierra de Juarez, Oaxaca, Mexico. Values are given for all species (dots, continuous line), and 
separately for terrestrials (downwards-pointing triangles, short-dashed lines) and epiphytes 
(upwards-pointing triangles, long-dashed lines). Trend lines were set with locally weighted 
regressions (LOWESS). 


and a Tukey post-hoc test (P< 0.001). We compared species richness per plot 
of our gradient with other similar transects from Bolivia, Ecuador, Costa Rica, 
and Mexico based on Kessler (2000, 2001), Kluge and Kessler (2005), Krémer et 
al. (2013), Salazar et al. (2013), and Carvajal-Hernandez et al. (2017). To test 
the differences in mean plot richness of all transects, we used an analysis of 
variance and a Tukey post-hoc test (P< 0.001). Because the general objective of 
this work and the general richness pattern was not altered, we included 
atypical data from the site at 2,500 m located on the other side of the mountain 
chain (leeward), which was excluded only in the comparison between 
transects. All data analyses were conducted in the statistical package R (R 
Core Team, 2016. 


RESULTS 


In the 121 plots, we recorded 195 species of ferns and lycophytes belonging 
to 60 genera and 21 families, including two new records for Oaxaca State, 
Ctenitis salvinii and Phlegmariurus capillaris, as well as one species recently 
described as new to science (Ceradenia sacksii; Sundue, 2017) (TABLE 1). The 
most species-rich genera were Elaphoglossum (20 species), Asplenium (18), 
and Pleopeltis and Hymenophyllum (12 each). The most species-rich families 
were Polypodiaceae (49), Dryopteridaceae (35) and Hymenophyllaceae (25). 

Fern species richness followed a hump-shaped elevational pattern with 
highest species numbers per plot at mid-elevations between 1,500 and 2,100 m 
and a marked decline towards both low and high elevations (Fig. 3). The 
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TaBLE 1. Species list of lycophytes and ferns recorded along the elevational gradient in the Sierra 
de Judrez, Oaxaca, Mexico. Data are presented with the minimum (Min) and maximum (Max) 
observed elevational distribution, frequency of the species (F; the number of plots in which the 
species was recorded), the mean number of individuals per plot (N), and substrate including the 
life form: t: terrestrial, he: hemiepiphytic, e: for epiphytes comprising the Johansson (1974) 
zonation in which the epiphytes occurred (from zone 1 to zone 5), dw: dead wood, and r: rocks. 
Species marked with " are new records for Oaxaca and with * were found outside or close to the 


investigated plots. 


Taxa 


Lycopodiaceae 

”Phlegmariurus capillaris (Sodiro) B. 
Ollg. 

*Phlegmariurus cuernavacensis (Underw. 
& F.E. Lloyd) B. Allg. 

Phlegmariurus dichotomus (Jacq.) 
W.H.Wagner 

Phlegmariurus linifolius (L.) B.Olig. 

Phlegmariurus myrsinites (Lam.) B. Ollg. 

Phlegmariurus pringlei (Underw. & F.E. 
Lloyd) B. Allg. 

Phlegmariurus taxifolius (Sw.) A. Live & 
D. Love 

Phlegmariurus wilsonii (Underw. & 
F.E.Lloyd) B.dllg. 

Selaginellaceae 

Selaginella aff. subrugosa Mickel & Beitel 

Selaginella flexuosa Spring 

Selaginella mosorongensis Hieron. 

Selaginella oaxacana Spring (GCI) 

Selaginella subrugosa Mickel & Beitel 

Selaginella sp. 

Aspleniaceae 

Asplenium abscissum Willd. 

Asplenium auriculatum Sw. 

Asplenium cirrhatum Rich. ex Willd. 

Asplenium cristatum Lam. 

Asplenium cuspidatum Lam. 

Asplenium fragrans Sw. ; 

Asplenium hallbergii Mickel & Beitel 

Asplenium harpeodes Kunze 

Asplenium miradorense Liebm. 

Asplenium monanthes L. 

Asplenium polyphyllum Bertol. 

Asplenium pteropus Kaulf. 

Asplenium salicifolium L. 

Asplenium radicans L. 

Asplenium serratum L. 

Asplenium uniseriale Raddi 

Asplenium sp.(cf. A. uniseriale Raddi) 

Hymenasplenium riparium (Liebm.) L. 
Regalado & Prada 


, Min 
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2800 
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1500 
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19 
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TABLE 1. 


Taxa 


Athyriaceae 

Diplazium diplazioides (Klotzsch & 
H.Karst.) Alston 

Diplazium franconis Liebm. 

Diplazium lonchophyllum Kunze 

Diplazium plantaginifolium (L.) Urb. 

Diplazium striatum (L.) C. Presl 

Diplazium ternatum Liebm. 

Blechnaceae 

Austroblechnum wardiae (Mickel & 
Beitel) Gasper & V.A.O.Dittrich 

Lomaridium ensiforme (Liebm.) Gasper & 
V.A.O.Dittrich 

Lomaridium fragile (Liebm.) Gasper & 
V.A.O. Dittrich 

Parablechnum falciforme (Liebm.) Gasper 
& Salino 

Cyatheaceae 

Cyathea aristata Domin 

Cyathea divergens Kunze 

Cyathea firma (Bak.) Domin 

Cyathea fulva (M. Martens & Galeotti) 
Fée 

Cyathea liebmannii Domin 

Cyathea myosuroides (Liebm.) Domin 

Cyathea salvinii (Hook.) Domin 

Cyathea schiedeana Domin 

Sphaeropteris horrida (Liebm.) R.M. 
Tryon 

Cyathea sp. 

Dennstaedtiaceae 

Histiopteris incisa (Thunb.) J. Sm. 

Hypolepis nigrescens Hook. 

Pteridium feei (W. Schaffn. ex Fée) Faull 

Dicksoniaceae 

Lophosoria quadripinnata (J.F.Gmel.) 
C.Chr. 

Dryopteridaceae 

Arachniodes denticulata (Sw.) Ching 

Bolbitis portoricensis (Spreng.) 
Hennipman 

Ctenitis sp. (cf. C. excelsa (Desv.) Proctor) 

Ctenitis hemsleyana (Baker ex Hems!l.) 
Copel. 

Ctenitis melanosticta (Kunze) Copel. 

Ctenitis nigrovenia (Christ) Copel. 

"Ctenitis salvinii (Baker) Stolze 

Dryopteris nubigena Maxon & C.V. 
Morton 

Dryopteris wallichiana (Spreng.) Hyl. 

Elaphoglossum affine (M.Martens & 
Galeotti) T. Moore 


Continued. 

Min Max. 
1000 2100 
2100 2100 
600 600 
600 600 
1000 1000 
1000 2100 
2100 2100 
1500 2200 
1500 2000 
1800 2200 
1500 1500 
1500 1500 
1500 1500 
1500 2200 
1000 1000 
1000 1000 
1500 2000 
1000 1500 
1000 1000 
2100 2100 
2000 2000 
1000 1000 
2500 2900 
2100 2200 
1300 2200 
500 500 
300 1300 
2000 2200 
500 500 
300 300 
600 600 
2800 3000 
2700 3000 
2800 2900 
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TABLE 1. Continued. 


Taxa Min Max. N F Life form 
Elaphoglossum alan-smithii Mickel 2100 2200 6.25 4 e1, e2, e3, 64, e5 
Elaphoglossum dissitifrons Mickel 2000 2900 11.71 17 t,e1, e2, e3, e4, ed 
Elaphoglossum glabellum J.Sm. 10007, 1000.91.00" 1" et 
Elaphoglossum glaucum T. Moore 1800 3000 13.41 22 4, 17).61ve2meRnes 
Elaphoglossum sp. (cf. E. glaucum T. 2100 2200 23.00. 2” e1;\e2; e3ne4 


Moore) 
Elaphoglossum guatemalense (Klotzsch) ‘1000 2000 7.43 21 t, dw, e1, e2, e3, e4, e5 
T. Moore 


Elaphoglossum leonardii Mickel 2800 2800 34.00 2 r 
Elaphoglossum lonchophyllum (Fée) T. 1500 2100 10.33, 03 e1,e2ned 
Moore 
Elaphoglossum muscosum (Sw.) T. 1500 1500) 54.00) ( wie? 
Moore 
Elaphoglossum paleaceum (Hook. & 2100 2200 6.64 11  e1;,62; e3;e4je5 
Grev.) Sledge 
Elaphoglossum parduei Mickel 2000 2000 5.80 5 e1,e2,e3 
Elaphoglossum peltatum (Sw.) Urb. 1500 2200 18.25 24 t, dw, rc, el, e2, e3, e4, e5 
Elaphoglossum petiolatum (Sw.) Urb. 2900 2900 1.00 2 e1, e2, e3, e4 
Elaphoglossum pringlei (Davenp.) C.Chr. 2500 2500 8:33 3° tdi 
Elaphoglossum sartorii (Liebm.) Mickel 600 2200 22.17 35 t, dw,r, el, e2, e3, e4 
Elaphoglossum seminudum Mickel 2200 2200 80.25: 4 ited . 
Elaphoglossum squamipes (Hook.) T. 2700 2900 4.00 3 (t), dw, r, e1, e2 
Moore 
Elaphoglossum tenuifolium (Liebm.) T. 2800 2900 7.80 5 tr 
Moore 
Elaphoglossum vestitum (Schltdl. & 600° 2100 17.57 37°61, 62)e3/¢4,/e5 
Cham.) Schott 
Elaphoglossum sp. 1000 1000 2.00 1 e2 
Lastreopsis effusa (Sw.) Tindale 600 .600- 2:00 iat 
Megalastrum mexicanum R.C. Moran & J. 600; 600 13.25 2% 
Prado 
Polystichum hartwegii (Klotzsch) Hieron 1800 2700 7.50 8 t 
Polystichum mickelii A.R. Sm. 1300; 1300 2:33) °3igeg 
Stigmatopteris longicaudata C. Chr. 1560; 1500\g63°334 iid 
Gleicheniaceae 
Diplopterygium bancroftii (Hook.) A.R. 1800; 1800 . 200%) 17% 
Sm. 
Hymenophyllaceae 
Abrodictyum rigidum (Sw.) Ebihara & 1000; 1800 26:738) 1ivet.= 
Dubuisson 
Didymoglossum hymenoides (Hedw.) 1300°- 1300 250), 2 Pei@ez 
Copel. 


Didymoglossum kraussii (Hook. & Grev.) 300. 500 ({.3:33:)) .9 4diwfetiad 
C. Pres] 


Didymoglossum ovale E. Fourn. 300 500 3/008) tearpe 

Didymoglossum reptans (Sw.) C. Pres] 1500 2100 10.00 12 t, dw,r, el, e2, e3 

Hymenophyllum asplenioides Sw. 1500 2200 4.88 8 (t),e1, e2 

* Hymenophyllum brevistipes Liebm. 1200) 1200 1:00°) 1 we2 

Hymenophyllum ectocarpon Fée 1800 3000 10.73 40 t, dw,r, e1, e2, e3 

Hymenophyllum fragile (Hedw.) C.V. 2100 2200 7.88 8 t,dw,r, e1, e2 
Morton 


ee 
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TABLE 1. Continued. 


Taxa Min Max. N F Life form 
Hymenophyllum fucoides Sw. 2100/'432200 1.441:88 68.8 4itpelye2 
Hymenophyllum hirsutum (L.) Sw. 1000 1500 3:33.01. 3itelziez,) a3 
Hymenophyllum Janatum Fée 15002200 4)-2:15 130) )(t}pel}e2;\e3 
Hymenophyllum myriocarpum Hook. 2000 2900 18.95 20 t, dw, e1, e2, e3, e4 
Hymenophyllum polyanthos (Sw.) Sw. 1000 2000 SSIS tpaw,y tel R62 793 


Hymenophyllum tegularis (Desv.) Proctor 2800 2800 4.00 1 r 
& Lourteig 


Hymenophyllum trapezoidale Liebm. 15002200 9:21. 14--: t,1r, 1, 62, e3 

Hymenophyllum tunbrigense Sm. 3000 3000 13.00 1 t,dw,ei 

Hymenophyllum undulatum (Sw.) Sw. 1500 1800 §.40/ 4 5), e1,.6e2;'68 

Polyphlebium capillaceum (L.) Ebihara& 1000 2100 2.80 5° dw, el, e2 
Dubuisson 


Polyphlebium hymenophylloides (Bosch) 1800 2200 2.86 7. t, dw, el, e2 
Ebihara & Dubuisson 

Polyphlebium pyxidiferum (L.) Ebihara & 500 1500 1.50 4 r,e1, e2,e4 
Dubuisson 


Trichomanes lucens Sw. 1800 2000 1.00 2 

Trichomanes polypodioides L. 10006 15007 143.00 |\"3 1,62 

Trichomanes crispum L. 1300 1300 S00; Lint red 

Vandenboschia collariata (Bosch) Ebihara 600 600 3.67 at ir, et 

& K. Iwats. 

Vandenboschia radicans (Sw.) Copel. 1300 1500 6.20 5 t, dw, r, el, e2 
- Lindsaeaceae 

Lindsaea arcuata Kunze 1000 1500 4:08.33 2 

Lindsaea lancea Christ 1000 1000 1.00 1s at 

Lomariopsidaceae 

Lomariopsis mexicana Holttum 1300 1300 2.00 1 t, he,(e1),( e2) 

Lomariopsis recurvata Fée 300 600 4.64 11 t, he, (r),(e1), (e2), (e3) 

Lygodiaceae 

Lygodium heterodoxum Kunze 300 300 1.50 2 t, he, (e1) 

Marattiaceae 

Danaea moritziana C. Pres] 1000 1300 7.50 Ae at 

Marattia excavata Underw. 1500 1500 1.00 Siar 

Marattia weinmanniifolia Liebm. 2100 2100 2:00 20 4Rt 

Plagiogyriaceae 

Plagiogyria pectinata (Liebm.) Lellinger 2800 3000 34.50 12 t, (e2) 

Polypodiaceae 

Alansmia cultrata (Willd.) Moguel & M. 1800 2900 19.64 14 r,e1, e2,e3 

Kessler 
Campyloneurum xalapense Fée 600 1500 10.22 18 t, dw, e1, e2, e3, e4, e5 


Campyloneurum amphostenon (Kunze ex 2100 2900 2.75 4 e2 
Klotzsch) Fée 
Campyloneurum angustifolium (Sw.) Fée 300 1000 5.73 22  t, e1, e2, e3, e4, e5 


Ceradenia oidiophora (Mickel & Beitel) 2100 2100 2.00 2° *e1, e2 
A.R. Sm. 

"Ceradenia sacksii Sundue 2800 2800 20.00 rs was 

Cochlidium linearifolium Maxon ex C. 1000 3000 37.11 37. e1, e2, e3, e4, e5 
Chr. 

Galactodenia delicatula (M.Martens & 2000 3000 5.77 22. e1, e2, 63, e4 


Galeotti) Sundue & Labiak 
Melpomene leptostoma (Fée) A.R. Sm. & 1500 2200 7.25 24. e1, e2, e3, e4, e5 
R.C. Moran 
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TABLE 1. Continued. 
Taxa Min Max N 

Melpomene moniliformis (Lag. ex Sw.) 2800 2900 63.73 
A.R. Sm. & R.C.Moran | 

Melpomene pilosissima (M.Martens & 2800 2800 5.00 
Galeotti) A.R.Sm. & R.C. Moran 

Melpomene xiphopteroides (Liebm.) A.R. 1000 2000 3.20 
Sm. & R.C.Moran 

*Moranopteris basiattenuata (Jenman) '2400 2400 1.00 
R.Y. Hirai & J.Prado 

Moranopteris trichomanoides (Sw.) R.Y. 1800 1800 2.00 
Hirai & J.Prado 

*Mycopteris semihirsuta (Klotzsch) 2400 2400 1.00 
Sundue 

Niphidium crassifolium (L.) Lellinger 300 1300 = 3.47 

Pecluma longepinnulata (E. Fourn.) F.C. 1300 2100 10.53 
Assis & Salino 

Pecluma atra (A.M.Evans) M.G. Price 300 300 Ae 

Pecluma consimilis (Mett.) M.G. Price 500 600 5.20 

Pecluma dispersa (A.M.Evans) M.G. Price 500 500 2.00 

Pecluma divaricata (E.Fourn.) Mickel & 500 500 ~~ 7.00 
Beitel 

Pecluma hartwegiana (Hook.) F.C. Assis 2700 3000 78.24 
& Salino 

Pecluma rhachipterygia (Liebm.) F.C. 600 1300 8.00 
Assis & Salino 

Pecluma sursumcurrens (Copel.) M.G. 500m7r 500 1.00 
Price 

Phlebodium areolatum (Humb. & Bonpl. 300)82200 © 3.16 
ex Willd.) J. Sm. 

Pleopeltis alansmithii (R.C. Moran) A.R. 2100 3000 33.25 
Sm. & Tejero 

Pleopeltis angusta Humb. & Bonpl. ex 600 1800 3.25 
Willd. 

Pleopeltis astrolepis (Liebm.) E. Fourn. 300 600 1.60 

Pleopeltis conzattii (Weath.) R.M. Tryon 2500 2500 1.00 
& A.F. Tryon | 

Pleopeltis crassinervata T. Moore 600 2200 7.18 

Pleopeltis fallax (Schltdl. & Cham.) 500" 96005, 11.67 
Mickel & Beitel 

Pleopeltis furfuracea (Schltdl. & Cham.) 500 (600 6:25 
A.R. Sm. & Tejero 

Pleopeltis madrensis (J.Sm.) A.R. Sm. & 2500 2500 £4.00 
Tejero 

Pleopeltis mexicana (Fée) Mickel & Beitel 1000 3000 34.68 

Pleopeltis muenchii (Christ) L.D. Gomez 2700 3000 56.15 

Pleopeltis plebeia (Schltdl. & Cham.) A.R. 1000 3000 18.06 
Sm. & Tejero 

Pleopeltis polypodioides (L.) E.G. 300 1500 8.00 
Andrews & Windham 

Pleopeltis pyrrholepis (Fée) A.R. Sm. & 300 1000 5.10 


Tejero 
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TABLE 1. Continued. 


Taxa Min Max. N Ke Life form 
Polypodium echinolepis Fée 1300 1500 5.18 11 e1, e2, e3, e4, e5 
Polypodium martensii Mett. | 2500 2900 34.86 21 (t), e1, e2, e3, e4, e5 
Polypodium plesiosorum Kunze 600 600 4.50 2 e2,e4,e5 


Polypodium pleurosorum Kunze ex Mett. 1800 2200 9.88 24 t,e1, e2, e3, e4, e5 
Polypodium puberulum Schltdl. & Cham. 2000 2200 5.69 16 e1, e2, e3, e4, e5 


Polypodium rhodopleuron Kunze 2100 2100 1.50 2 t,e1, e2, e3, e4, e5 
Polypodium subpetiolatum Hook. 2700 2700 4.00 2 e1 
Polypodium sp. 1000 1000 3.00 1 e1 
Serpocaulon falcaria (Kunze) A.R. Sm. 1500 32700 | Zaa0 °28 t, Tel, eZ, :e3, 64 
Serpocaulon triseriale (Sw.) A.R. Sm. 300 1000 3.63 8 tt, dw, el, e2, e3, e4, e5 
Stenogrammitis delitescens (Maxon) 2200" 2200 F £6,00** 1° e2 
Labiak 
Stenogrammitis hellwigii (Mickel & 2800 3000 21.13 15 r,e1, e2 
Beitel) Labiak 
Stenogrammitis prionodes (Mickel & 1500 2200 9.33 18 e1, e2, e3, e4, e5 
Beitel) Labiak 
*Terpsichore aspleniifolia (L.) A.R.Sm. 2400 2400 1 Le; 
Pteridaceae 
Adiantum poiretii Wikstr. 2500, 2500. 7,00». .2),t,e1 
Adiantum pulverulentum L. S00 ee G00 pee 35- 0758 
Adiantum tenerum Sw. DOU. OU ere OD a kes t 
Adiantum trapeziforme L. SO0s -  OMUsF wct,00. a. Lat 
- Ananthacorus angustifolius (Sw.) 500! 500+) 43.50’. 2/e3 
Underw. & Maxon 
Antrophyum ensiforme Hk.in Benth. 1000 2100 6.36 11 e1, e2, e3, e4 
*Gaga chaerophylla (M. Martens & 1800 1800 1.00 1 t 
Galeotti) F. W. Li & Windham 
Jamesonia hirta (Kunth) Christenh 20005-2000 1,00 .1 % 
Polytaenium chlorosporum (Mickel & 500 ©1800. \42,.40 >" Stel 
Beitel) E.H.Crane 
Polytaenium feei (W.Schaffn. ex Fée) 500 600 11.40 5 (t), dw,r, el, e2, e3 
Maxon 
Pteris altissima Poir. 300 600 3.00 4 t,dw 
Vittaria bradeorum Rosenst. 1500 2200 2.82 11 dw, el, e2 
Vittaria flavicosta Mickel & Beitel 500 1000 2.89 9 (t) e1, e2, e3, e4 
Vittaria graminifolia Kaulf. 2100 3000 30.87 39 (t) dw,r, e1, e2, e3, e4, e5 
Saccolomataceae 
Saccoloma inaequale (Kunze) Mett. AOOO .e 2100 6 127 1a 2h et 
Tectariaceae 
Tectaria heracleifolia (Willd.) Underw. 300 600 50.63 19 t,r 
Thelypteridaceae 
Goniopteris biolleyi (Christ) Pic. Serm. BOOST GOO! LOOM Mitt 
Goniopteris ghiesbreghtii (Linden) J. Sm. GOO: 61 600s) 600 up <2 oth 
Goniopteris hatchii (A.R.Sm.) A. Léve & 1500 1500 12.50 2 t,dw 
D. Love 
Goniopteris imbricata (Liebm.) A. Love & 500 ~ 500”” 3.67... 6 “t 
D. Love 
Goniopteris liebmannii (Maxon & C.V. 200K SOCMINZI a FEA 


Morton) Salino & T.E. Almeida 
Goniopteris obliterata C. Pres] S00. OOM 65.676. situs 
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highest total species richness per elevational step was 57 species at 1,500— 
2,100 m and the lowest at 3,000 m with 15 species. 

Within the total of 195 species we found 141 growing as terrestrials and 129 
as epiphytes, whereas 75 were strictly terrestrial and 77 strictly epiphytic, 
which means that 39 species were recorded within both life forms, excluding 
four hemiepiphytic species. Species with double notations (as terrestrials and 
epiphytes) mainly belonged to the families Dryopteridaceae, Polypodiaceae, 
and Hymenophyllaceae. Both major life forms followed different trends, with 
epiphytic richness per plot mostly about double of the terrestrials, and with a 
pronounced hump-shaped peak, whereas terrestrials were more.or less 
constant up to 2,500 m and decreasing slightly at higher elevations (Fig. 3). 

Epiphytic species on average had a wider elevational range and frequency 
than terrestrials (elevational range of epiphytes: 510.3 + 580.8 m and 
terrestrials: 222.8 + 351.9 m; Wilcoxon test, P<0.01, frequency of epiphytes: 
9.61 + 10.06, both: 10.35 + 10.21 and terrestrials: 5.77 + 6.85, ANOVA, F= 
4.8, df =2, P < 0.01). 

All other transects from Bolivia, Ecuador, Costa Rica, and Mexico showed a 
hump-shaped richness pattern and similar maximum richness at mid- 
elevations (more or less at the same elevation). Only the transects located at 
the edges of the tropical realm, in Mexico and the southernmost transect, at 
Masicuri in Bolivia, presented lower richness values compared with Costa 
Rica, Ecuador, and the other transect in Bolivia (Fig. 4, ANOVA, F=20.27, df= 
7, P < 0.001). However, the peak of maximum richness did not decrease with 
increasing latitude towards lower elevations. , 


DISCUSSION 


With this study, we report the most detailed plot-based study on fern species 
richness along a complete elevational gradient (300-3,000 m) in humid 
tropical forest in Mexico. In 121 plots we found 195 species mostly belonging 
to the families Polypodiaceae, Dryopteridaceae, and Hymenophyllaceae, 
coinciding with other transects in Mexico (Acebey, Krémer, and Kessler, 
2017; Carvajal- Hernandez et al., 2017; Krémer et al., 2013), Central and South 
America (Kluge and Kessler, 2005; Mehltreter, 1995; Salazar et al., 2013). The 
species recorded in this study represent 29% of the total ferns of the state of 
Oaxaca in a surface of ca. 95,360 km? and 18% of Mexico (Mickel and Smith, 
2004), this in only 4.8 ha of sampling. This shows the importance of the Sierra 
de Juarez, as a national refuge for ferns, since in this sampling area occur more 
species than in larger areas as the Yucatan Peninsula (ca. 125,000 km’; 
integrated for three complete states in Mexico) with only 65 species 
(Fernandez Carnevali et al., 2012) or the state of Tabasco (ca. 24,731 km’) 
with 169 species (Carvajal-Hernandez et al., 2018). 

The highest species numbers were recorded at mid-elevations at about 
1,500—2,100 m. This hump-shaped pattern of fern richness is commonly found 
in tropical mountains and has been explained by climatic factors, especially 
temperature, precipitation, and fog frequency, e.g., from the Andes (Kessler, 
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Fic. 4. Left: Patterns of fern species richness along eight elevational transects in the Neotropics, 
following the latitudinal gradient from Bolivia (-19 and-16° S) via Ecuador (~0°) and Costa Rica (9° 
N) to Mexico: Los Tuxtlas (17° N) and Cofre de Perote (19° N), Veracruz and Sierra de Juarez, 
Oaxaca (18° N). Trendlines were set with locally weighted regressions (LOWESS). Right: Mean plot 
richness plotted against latitude. Letters above indicated significant different groups following 
analysis of variance (ANOVA) and Tukey post-hoc test (P< 0.001). 


2001; Krémer et al., 2005; Salazar et al., 2013), Central America (Watkins et al., 
2006), Asia (Bhattarai, Vetaas, and Grytnes, 2004; Grytnes and Beaman, 2006), 
and Africa (Hemp, 2002). However, while species richness counts, especially 
in the Neotropics, average 40-50 species in plots of the same size as in this 
study, we found a considerably lower number of fern species at mid-elevations 
(on average 25 species), which was comparable to nearby transects in Veracruz 
(Acebey et al., 2017; Carvajal-Hernandez et al., 2017; Krémer et al., 2013). 
Higher species richness is expected in riparian and ravine zones in the Sierra 
de Juarez, where ferns are more abundant due to humid conditions (Carvajal- 
Hernandez and Krémer, 2015; Kessler, 2010; Kluge and Kessler, 2011). 
Because this study excluded these habitats, the total number of species is 
likely to be somewhat higher than found here. 

In evaluating life forms, we observed that the elevational hump of species 
richness was mainly driven by epiphytes. Apart from the extremes at the upper 
and lower ends of the gradient, this life form had twice the number of species 
per plot as terrestrial taxa. This has also been reported before (Acebey, Krémer 
and Kessler, 2017; Kessler, 2001; Kluge and Kessler, 2005; Krémer et al., 2013; 
Salazar et al., 2013; Watkins et al., 2006) and seems to be a consistent tropics- 
wide fern richness pattern. In our study, we additionally found that epiphytes 
had a wider elevational range and were more frequent compared to terrestrial 
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species, as also found for ferns in Bolivia (Kessler, 2002). Dispersal might play 
a role in the broader distribution of epiphytes, because spores released in the 
forest canopy are likely to disperse greater distances (Peck, Peck, and Farrar, 
1990). Additionally, epiphytic and terrestrial niches operate at largely 
different scales: whereas terrestrial microhabitats change due to, e.g., soil 
nutrients, soil moisture, and topographic parameters over distances of several 
dozen meters (Jones, Szyska, and Kessler, 2011; Cicuzza et al., 2013), epiphytic 
niche differentiation is driven by small-scale features due to tree crown 
architecture (Kromer, Kessler, and Gradstein, 2007). Hence, a study plot of 400 
m* is likely to contain a wider range of epiphytic than terrestrial microhabitats 
with the consequence that this plot usually contains a higher number of 
epiphytes with a wider ecological variability. 

Probably the most interesting result of our study is that, at least for 
neotropical ferns, there appears to be no gradual downward shift in species 
richness with increasing distance from the Equator, but a sudden reduction of 
species richness at the transitions from the tropics to the subtropics, as already 
reported by Salazar et al. (2013). The causes of this sudden reduction in fern 
richness are probably related to climatic conditions but may also include a 
historical and evolutionary signal related to the tropical origin of many major 
fern groups (Kessler, Karger, and Kluge, 2016). 

Although the regional fern diversity of humid montane forests in 
southeastern Mexico is high, we found that at the local scale it is lower 
compared with similar ecosystems in Costa Rica, Ecuador, and Bolivia. The 
confirmation of the Mexican Biogeographic Transition Zone (MTZ) for ferns 
suggests that this zone is likely a consequence of climatic factors. Colder 
temperatures and especially the occurrence of frost act as climatically 
limiting factors of fern distributions with increasing elevation (Bhattarai, 
Vetaas, and Grytnes, 2004; Kluge, Kessler, and Dunn, 2006) and may be 
limiting factors in the latitudinal distributions of ferns with increasing 
distance from the Equator. Moreover, the lack of appropriate data from the 
latitudinal range between Costa Rica and Mexico indicates that the change in 
species richness may also either be abrupt further south or even more gradual 
across Central America. Therefore, more research is needed inside of this 
transition zone and in Central America to have a complete overview of the 
latitudinal diversity patterns. 
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APPENDIX 1. 
Mexico. 


Locality, 
Municipality 


Cerro Tejon, Paso Nuevo la 
Hamaca, Valle Nacional 
Cerro Marin, Valle Nacional 

Cerro Quemado, Vega del 
Sol, Jacatepec 

Yetla, Valle Nacional 

Puerto Antonio, Comaltepec 


Cerro Redondo, La 
Esperanza, Comaltepec 
Santa Cruz Tepetotutla, San 

Felipe Usila 
San Isidro, San Pedro Yolox 


Santa Cruz Tepetotutla, San 
Felipe Usila 

Jatianguis, Ixtlan 

Cerro Zacate, San Pedro 
Yolox 

Cerro Pel6n, Comaltepec 

Cerro Zacate, San Pedro 
Yolox 

Unknown, San Pedro Yolox 


Elevational 
band (m) 


300 


500 
600 


1,000 
1,300 


1,500 
1,800 
2,000 
2,100 
2,200 
2,500 


2,700 


2,800 
2,900 


3,000 


Coordinates 


17.78685 N, 96.23103 W 


17.82863 N, 96.25282 W 


17.79322 N, 96.17304 W 


17.72213 N, 96.3089 W 
17.67064 N, 96.33496 W 


17.62129 N, 96.34547 W 
17.71802 N, 96.5587 W 

17.60485 N, 96.41625 W 
17.68011 N, 96.55075 W 
17.68041 N, 96.55442 W 
17.3804 N, 96.51286 W 

17.61111 N, 96.53198 W 


17.58451 N, 96.51245 W 
17.61381 N, 96.53339 W 


17.60825 N, 96.52305 W 
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Location of the study sites along the elevational gradient in Sierra de Juarez, Oaxaca, 


Type 
of forest 


Rain Forest 


Rain Forest 
Rain Forest 


Rain Forest 

Transition Rain Forest — 
Cloud Forest 

Cloud Forest 


Cloud Forest with 
Oreomunnea mexicana 

Monodominant 
Oreomunnea Forest 

Cloud Forest 

Cloud Forest 

Pine Forest 

Oak Forest 


Oak forest 
Pine-Oak Forest 


Elfin Forest (tree line) 
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AsstTract.—Ferns form a part of vegetable consumption in many communities across the globe. 
Their popularity is rising because of their high nutritional content such as vitamin C, protein, and 
iron. Given that the nutritive value of different wild edible ferns has not been explored 
comprehensively, the present study was carried out. Six different species of wild edible ferns viz. 
Diplazium esculentum, Dryopteris cochleata, Diplazium spectabile, Diplazium dilatatum, Tectaria 
coadunata, and Diplazium maximum were analyzed for dietary nutrients and multi-elemental 
profiling. The highest moisture content (89.83%) and highest vitamin C content (25.38 mg g '°”) 
was found in Diplazium maximum.The highest protein was found in Tectaria coadunata (12.10%), 
which was much higher compared to many common commercial vegetables. Fat content was 
highest in Dryopteris cochleata (0.40%), whereas fiber content was found to be maximum in 
Diplazium esculentum (5.05%) and ash content was highest in Diplazium spectabile (1.57%). The 
mineral analysis of these ferns showed that they are good suppliers of sodium (0.46— 3.56 mg g '°°), 
magnesium (6.42— 13.39 mg g '°°), and iron (10.36— 16.35 mg g '°°). Most of the toxic and heavy 
metals were found to be within permissible limits. However, some of the fern species collected 
from certain regions had higher levels of these metals. Based on the present study on dietary 
nutrients and elemental profile, these edible ferns could be popularized for commercial 
cultivation. Further, research work could be conducted for standardizing practices and, because 
of anticipated climate changes, these ferns could add diversity in the world food basket. 


Key Worps.—elemental profile, green leafy vegetables, heavy metals, ICPMS, wild vegetables. 


Among all vegetables, green leafy vegetables have occupied a unique place 
because of their color, flavor, and health benefits. They are rich sources of B- 
carotene, ascorbic acid, iron, zinc, folate, and dietary fiber (Negi and Roy, 
2000). The nutritional values and mineral content of wild edible plants were 
richer than that of commercial vegetables and could be used for nutritional 
purposes (Seal, 2012). Wild plants serve as an indispensable constituent of the 
human diet, supplying the body with minerals, vitamins, and certain hormone 
precursors, in addition to protein and energy (Akubugwo, Obasi, and Ginika, 
2007). These plants also provide some useful products e.g., medicine, fiber, 
and fodder (Kayang, 2007). The study of these underutilized vegetables is 
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important not only to identify the potential sources that could serve as 
alternative food but also to select promising types for domestication. In recent 
years various lesser known plants are being evaluated for their nutritional and 
medicinal properties (Bassey et al., 2001; Caldwell, 2010; Pradhan et al,. 2015) 
but there are still large groups of species that are commonly eaten and sold at 
local market and whose nutritional potential has not been studied. Among 
wild plants, edible ferns play an important role as they are often included in 
the daily diet by the population of the Eastern Himalayas. 

Due to its varied topography, climate, soil, and vegetation, the Himalayan 
region is very rich in ferns and lycophytes. Out of 362 ferns belonging to 108 
genera available in the Sikkim Himalayas (Kholia, 2010) croziers of six species 
viz. Dryopteris cochleata (Buch.-Ham. ex D. Don) C. Chr., Diplazium spectabile 
(Wall. ex Mett.) Ching, Diplazium maximum (D. Don) C. Chr., Diplazium 
dilatatum Blume, Diplazium esculentum (Retz.) Sw., and Tectaria coadunata 
(J. Sm.) C. Chr. are popularly consumed as vegetables. The versatility of their 
taste makes them a good selection for vegetable purpose. In the Sikkim 
Himalayas, croziers of these six species are sold in every local market during 
each fern’s peak season, which varies from March to December. These ferns are 
considered to possess high amounts of dietary nutrients and other secondary 
metabolites. However, less effort has been made towards the research and 
development of these plant species. 

Even when research was conducted, it focused primarily on D. esculentum 
(Irawan et al., 2006; Saha, Biswal, and Deka, 2015; Seal, 2012; Tag et al., 2014; 
Tongco et al., 2014). Therefore, the present research was carried out to 
discover the nutrient potential of other wild leafy edible ferns. 


MATERIALS AND METHODS 


Diplazium esculentum, D. cochleata, D. spectabile, D. dilatatum, T. 
coadunata, and D. maximum were the experimental materials in the present 
study. The research was carried out during the years 2014-2016 in the P.G 
Laboratory, Department of Horticulture, Sikkim University, Gangtok, India. 

Collection of samples.—Fresh young samples (unopened young ferns leaves 
[croziers] of 20 to 25 cm length) were collected from different localities in the 
Sikkim Himalayas (Fig. 1) using GPS (Garmin) during the months of June— 
November 2015 and February—April 2016 (Table 1). The ferns were then stored 
in fresh sterile polyethelene bags at —20°C until analysis. The identification of 
the specimens was confirmed by consulting available literature published on 
ferns and through expert service from the Botanical Survey of India (BSI, 
India). 

Determination of dietary nutrients.—Moisture content of the leaves was 
estimated by the gravimetric method. Crude fat was determined by essential oil 
extractor, (Socsplus-SCS 06 DLS, PELICAN, India) by following the solvent 
extraction method (AOAC, 2000). Crude fiber was analyzed using the fiber 
estimation system, (Fibra plus-FES 04 AS DLS, PELICAN, India). In addition, 
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MR 1540s 1610 ° 
YNII 1627 KH 1132 "UN 1220 
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AR 11522 


LI 2350 


AL- Assam Lingey; AL- Assam lingey; AR |- Aritar; AR II- Aritar Il; BH- Bhusuk; 
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KA- Kabi; KH-Khamdong; L |- Lava | (kalimpong); LII- Lava 1; LN- Lower Nandok; 

el ee a a MR- Martam; NA- Nandok; SM- 6th mile; SO- Sombaria; THM- 32 Mile; UN- Upper Nandok; 
YN I- Yangyangl; YN II-Yangyang Il; YNIII- Yangyang Ill. 


Fic. 1. Location details of fern collection sites. 


ash content was estimated by placing a crucible with dried sample in an 
electrical muffle furnace (Yorco, India) at 600°C for 6 hours. 

Crude protein was determined by Lowry’s method (Thimmaiah, 2004) using 
a UV/VIS Spectrophotometer (Perkin Elmer, Lambda 35 UV/VIS spectrometer). 
Vitamin C was determined by titremetry (Rangana, 2012). 

Analysis of ionomic profiling was carried out by the ICPMS (Perkin-Elmer 
Nex ion 300X) with cross nebulizer. Digested samples were analyzed for the 
ionic constitution after standardization of the instrument using the multi 
elemental standard solutions supplied by the manufacturer for Ag, Al, As, Ba, 
Be, Bi, Ca, Cd, Ce, Co, Cs, Cu, Fe, Ga, Hg, I, In, Li, Mg, Mn, Mo, Na, Ni, Pb, Rb, 
Si, Sr, Ti, U, V, Xe, and Zn analytes. 

Statistical analysis.—The treatments were subjected to ANOVA using 
Completely Randomized Design (CRD) for dietary nutrient analysis using the 
software STPR. Each species was considered as one treatment and the different 
collection localities for each species were considered as replications. The 
number of treatments was six and the replications were four each per 
treatment. 


RESULTS 


Ferns were analyzed for proximate constituents viz. moisture content, crude 
fat, fiber, ash, and vitamin C as well as for protein, and the results (Table 2) 
revealed that edible ferns could act as a supplement for many nutrients in the 
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TABLE 1. Location details of Accessions. 
Species with Altitude 
harvesting time Places Latitude Longitude (m MSL) 
Dryopteris cochleata Yangyang III N 27°18.244’ E088°24.996' 1467 
(March-June) Aritar NI27°11,176" E088°40.552’ 1522 
Lower Nandok N 277177031, E088°36. 443’ 1174 
Assam Lingey I IN 27-7 E088°37.30' 966 
Tectaria coadunata Nandok N 2770715: E088°36. 837’ 1220 
(April - July) Dochum N 27°14.508' E088°28.942’ 366 
Sombaria N 27°10.993’ E088°11.864' 366 
Bermoik N 27°14.608' E088°28.337 531 
Diplazium dilatatum Bhusuk N 27°35.087’ E088°67.658' 1132 
(August-November) Dzongu N 27°27.40' E088°30.42' 1012 
Lava I N 27°12085’ E088°56.09' 2350 
Kabi N 27°24.501’ E088°37.561' 1443 
Diplazium esculentum Bermoik N 27°14.608' E088°28.337' 531 
(March- June) Khamdong N 27°17.68' E088°30.42’ 1132 
Yangyang I N 27°27.68' E088°25.18' 1303 
Bhalukhop, kpg N 27°05.522’ E088°29.419’ 997 
Diplazium spectabile Martam N 27°29.035' E088°55.196' 1540 
(August- November) 32 Mile, Gangtok  N 27°15.822’ E088°33.749' 617 
Upper Nandok N&2y 17715" E088°36. 837’ 1220 
6th mile N 27°18.495/ E088°035.307' 1610 
Diplazium maximum Yangyang II N 27° 17.661’ E088°24.789' 1527 
(October - December) Lava II N 27°107.741' E088°56.096' 2350 
Assam lingey II N 27°15.35’ E088°30.653' 966 
Aritar N 27°11.245/ E088°40.428' 1500 
TABLE 2. Proximate nutrient content of wild edible fern and its allies. 
Moisture Crude Fat Fiber Ash Vitamin C Protein 
Species (%) (%) (%) (%) (mg g *°°) (%) 
Diplazium 87.38 (69.22). 0.18 (2.38) 4.93 (12.81) 1.46 (6.86) 21.25 6.08 (14.05) 
dilatatum 
Diplazium 89.34 (70.90) 0.25 (2.74) 5.05 (12.98) 1.33 (6.61) 21.38 3.84 (11.11) 
esculentum 
Diplazium 89.83 (71.40) 0.20 (2.52) 3.91 (11.33) 0.98 (5.65) 25.38 5.58 (13.42) 
maximum 
Diplazium 88.16 (69.83) 0.21 (2.61) 3.39 (10.57) 1.57 (7.10) 21.63 3.56 (10.65) 
spectabile 
Dryopteris 82.23 (66.19)* 0.40 (3.64) 3.77 (11.16) 1.10 (5.99) 23.13 11.84 (20.10) 
cochleata 
Tectaria 85.51 (67.39) 0.18 (2.38) 2.35 (8.74) 1.23 (6.33) 22.25 12.10 (20.27) 
coadunata 
GM 87.07 (69.15) 0.24 (2.71) 3.90 (11.26) 1.28 (6.42) 22.50 7.16 (14.93) 
CD at 5% 1.60 (1.57) 0.13 (0.75) 0.63 (0.98) 0.22 (0.60) 2.26 2.48 (3.14) 


* (The value in parenthesis is in arcsine square root transformed values) 
(GM- Grand mean) 
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Tass 3. Concentration of dietary minerals in wild edible ferns (mg g '°°) (Mean with SD values). 


Species Ca Cu Fe I Mg Mn Mo Na Zn 
Diplazium aides G19, 2055. 610,622... 0,232),, 13.392). ,.6,542 J 96S hl O11S 
dilatatum 0.14 Foodies 510 0.03 6.28 4.00 7.70 1.27 0.03 
Diplazium mae 13.07 2) 4 So 210,195 9.065 ites 0.192 0.502 0.10 
esculentum 0.22 5.07 16.88 0.09 3.83 14.08 0.09 0.10 0.01 
Diplazium apemerig7os 16.352 50/172 6.42+ Liu U9" 3.562 2.045 
maximum 0.73 17.09 17.00 0.05 4.24 5.08 1.25 6.16 5.42 
Diplazium Peeve 25102 810-3626 0,23 2) °12,002 s1isi2t  ».0.482,..0.46%) 0.23+ 
spectabile 0.23 2.14 6.57 0.03 2.86 5.39 0.37 0.08 0.17 
Dryopteris Peewee cece 10.902 4/0192 “12.9/5 20.35= 0.10: (0.502 0.092 
cochleata AD a #- 3.03 7.50 0.05 4.76 6.16 0.02 0.12 0.03 
Tectaria Peerage rh 0lS 10255 Sales = 12.002" 0.162 ~ 6.55%°" .0,12> 


coadunata 0.18 2.83 3:85 0.01 3.63 3.54 0.08 Oftd 0.03 


human diet. Other than D. esculentum, this is the first report on proximate 
components and other dietary nutrients of wild edible ferns of eastern 
Himalayas. 

Moisture content.—The results showed that the fronds had high moisture 
content ranging from 82.23% (D. cochleata)—89.83% (D. maximum). 

Crude fat.—The fat content analyzed ranged from 0.18% (T. coadunata and 
D. dilatatum) —0.40% (Dryopteris cochleata) and Dryopteris cochleata was 
significantly superior to other species. 

Crude fiber content.—Diplazium esculentum was found to have superior 
amounts (5.05%) of crude fiber content followed by D. dilatatum (4.93%). 
These two species had significantly superior fiber content to the rest of the 
species. T. coadunata has the smallest amount (2.35%) of fiber. 

Ash content.—Ash content ranged from 0.98% (D. maximum) to 1.57% (D. 
spectabile). Diplazium esculentum had ash content of 1.33%. 

Vitamin C.—The highest amount (25.38 mg g '°°) of vitamin C was found in 
D. maximum, whereas D. cochleata (23.13 mg g '°°) was on par with D. 
maximum. The smallest amount (21.25 mg g '°°) was found in D. dilatatum. 

Protein.—Significantly higher amounts of protein content were found in T. 
coadunata (12.10%) and D. cochleata (11.84%) compared to other species. The 
smallest amount (3.56%) was found in D. spectabile. 

Dietary minerals.—Multielemental profiling was carried out for all the fern 
species collected. We analyzed 32 different elements. Only elements that have 
dietary importance (Table 3) and toxicity to humans (Table 4) are discussed 
here. National Institute of Nutrition Research, Hyderabad, India has listed 12 
elements for human dietary needs (Gopalan, Sastri, and Balasubramanian, 
1989) of which three elements were not analyzed and the remaining nine 
elements are reported here. 

The calcium content in the present experiment revealed that it ranged 
between 0.66 mg g ‘°° in D. esculentum—-1.09 mg g ‘°° in T. coadunata. The 
copper content was estimated to be lowest in T. coadunata (11.29 mg g '°”) 
and highest in D. dilatatum (19.25 mg g '°°). The present study revealed the 
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Tas_e 4. Concentration of heavy metals and toxic elements in wild edible ferns (ppm) (mean with 
SD values). 


Species Al As Cd Hg Li Ni Pb 
Diplazium dilatatum 52:9 ues Gost 0.01+ He 20,783 O2+ 
24.4 3.4 0.2 0.01 0.2 13.3 0.1 
Diplazium esculentum 58.5 14.6+ 0.4+ 0.07 = 24 24.5+ 06+ 
29.2 7.6 0.2 0.1 uo 10.8 0.6 
Diplazium maximum 291575 1633 ~ Ate 0.4+ 25.5 29.5> D+ 
324.1 7.4 14.4 0.7 29.4 20.1 0.1 
Diplazium spectabile 89152 16.2 0.6+ 0.01+ 11s 26+ i+ 
oe 8.7 0.2 0.02 0.2 oN 2:1 
Dryopteris cochleata 116.9+ 18+ 0.4+ 0.02+ 1.64 57.54 a3+ 
53.4 0.01 0.1 0.01 0.4 Ps | 0.2 
Tectaria coadunata 129.3% 1 SS fe 2.02 01> 193 32.84 0.5 
29.5 10 5.4 0.01 0.7 18.2 0.4 


iron content in most of the fern species analyzed was usually in the high range, 
from O:36sme soe. ld spectabile) -16.35 mg prey (D. maximum). Three 
species namely D. dilatatum, T. coadunata, and D. spectabile had the highest 
amount (0.23 mg g '°°) of iodine whereas D. maximum was estimated to have 
the smallest amount (0.17 mg g ‘’’). In the present study, magnesium ranged 
from 6.42 mg g ‘°° (D. maximum) -13.39 mg g'‘” (D. dilatatum). The 
manganese content ranged from 6.54 mg g-'”° (D. dilatatum.) — 20.35 mg g *°° 
(D. cochleata). Molybdenum content was estimated to be highest in D. 
dilatatum (3.96 mg g ‘°°) followed by D. maximum (0.99 mg g*%) and D. 
spectabile. (0.48 mg g ‘°°) and lowest in D. cochleata (0.10 mg g '°°) whereas, 
sodium content was found to be maximum in D. maximum (3.56 mg g '"”) and 
minimum in D. spectabile (0.46 mg g-'°°). Among the six edible ferns D. 
cochleata (0.09 mg g ‘°°) had the lowest value for zinc and D. maximum (2.84 
mg g -”) had the highest value. 

Toxic and heavy metals —The aluminium content in the edible ferns 
showed relatively higher values. Diplazium maximum was estimated to 
possess maximum Al content (291.7 ppm) followed by T. coadunata (129.3 
ppm).The arsenic content varied from 14.6 ppm (D. esculentum)—18.7 ppm (D. 
dilatatum) in the analyzed species. The cadmium content was highest in D. 
maximum (7.9 ppm) and lowest in D. esculentum (0.4 ppm). The mercury 
content ranged between 0.01 ppm (D. dilatatum, D. spectabile, and T. 
coadunata)— 0.07 ppm (D. esculentum). Lithium content ranged from D. 
maximum (25.5 ppm)—1.1 ppm (D. spectabile), while lead content varied from 
0.2 ppm (D. dilatatum) -1.3 ppm (D. spectabile). Nickel content was in the 
range of 20.7 ppm (D. dilatatum)-57.5 ppm (D. cochleata). 


DISCUSSION 


Moisture content.—Diplazium esculentum had a moisture content of 
89.34%, which was similar to the findings of earlier workers with 90.84% 
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(Irawan et al., 2006) and 91.82% (Tongco et al., 2014). Since there has been no 
work on other fern species, they were compared with commonly available 
leafy vegetables (Table 5). The moisture content of other leafy vegetables 
ranged from 78.66% (Drumstick)—95.43% (Celery). In the present study the 
moisture content of ferns varied from 82.23% (D. cochleata)—89.83% 
(Diplazium maximum), which was comparable with most of the leafy 
vegetables. The range was nearer to taro leaves (85.66%) and tree fern 
(88.60%). 

Fat content.—Fat content in the present study was in the range of 0.18%-— 
0.40% which were comparable to Amaranthus (0.33%), asparagus (0.12%), 
celery (0.17), coriander (0.52%), lettuce (0.15), ostrich fern (0.40%), and 
spinach (0.39%). Only taro leaves (0.74%) and drumstick leaves (1.4%) had 
higher fat content than the analyzed ferns. Earlier reports on the fat content of 
the D. esculentum had reported it as 0.34% (Saha, Biswal, and Deka, 2015) and 
0.36% (Pradeepkumar et al., 2015). However, according to Handique (1993) 
young D. esculentum contained very low fat content (0.18%). In the present 
study, the fat content of D. esculentum was estimated at 0.25 %. 

Fiber content.—Diplazium esculentum had highest fiber content (5.05%). 
Earlier research on D. esculentum by Irawan et al., (2006) and Tongco et al., 
(2014) reported crude fiber to be 4.82%. In addition Saha, Biswal, and Deka 
(2015) also reported in the same species the crude fiber to be 4.45%. Among 
the many leafy vegetables used in the human diet, only Amaranthus viridis 
(6.1%) and curry leaves (6.4%) (Gopalan, Sastri, and Balasubramanian, 1989) 
had higher fiber content than the D. esculentum. All other fern species also had 
a good amount of fiber, with the lowest being in T. coadunata (2.35%). All 
other leafy vegetables had fiber content from 0.50% (water cress) to 3.70% (taro 
and tree fern). Hence, all the ferns analyzed had higher fiber content. 

Ash content.—Based on the range of ash content (0.98-1.57%) in the ferns 
that were analyzed, they were as good as Amaranthus (1.5%), coriander 
(1.47%), and water cress (1.20%). In contrast, asparagus, celery, lettuce, and 
tree ferns had the smallest amount of ash content. However, drumstick, 
spinach, and taro leaves had higher content among leafy vegetables than the 
ferns. 

Vitamin C.—A comparison with some of the leafy vegetables revealed that 
the vitamin C content in ferns could be comparable to spinach (28 mg g '°°) 
and coriander leaves (27 mg g ‘°’). Many of the leafy vegetables viz. asparagus 
(5.6 mg g °°), celery (3.1 mg g '°°), and lettuce (9.2 mg g °°) had lower 
vitamin C content than the ferns. The other edible ferns such as Ostrich fern 
(26 mg g '°°) and tree fern (30 mg g '"’) had similar values to the ferns 
surveyed in the present study. Vitamin C content in D. esculentum (21.38 mg 
g '°°) was higher than the earlier finding of Irawan et al., (2006). But other 
reports estimated it to be even higher at 29 mg g ‘°° (Handique 1993) and 32.03 
mg g ‘°° (Saha, Biswal, and Deka, 2015). 

Protein content.—Ferns in the current study contained high amounts of 
protein. The highest content was estimated in Tectaria coadunata (12.10%) 
and Dryopteris cochleata (11.84%). These values were higher than the protein 
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content of drumstick leaves (9.40%), taro leaves (4.98%), Amaranthus (2.46%) 
or any of the other leafy vegetables under comparison. Even D. spectabile, 
which had the lowest amount (3.56%) of protein, was better than many leafy 
_ vegetables such as amaranth, asparagus, celery, lettuce, and spinach. 
Comparison with other earlier reported edible ferns revealed that the 
Himalayan ferns, except D. spectabile and D. esculentum, generally had a 
higher amount of protein than the ostrich fern (4.55%) and Marsilea crenata 
(Amoroso, Mendez, and Villalobos, 2017). The protein content in T. coadunata 
and D. cochleata was higher than most of the common vegetables such as 
cabbage (1.8%), cauliflower (5.9%) and even leguminous vegetables such as 
peas (7.2%) and broad bean leaves (5.6%) (Gopalan, Sastri, and Balasubrama- 
nian, 1989). These commercial vegetables are available during the winter, but 
the edible ferns can be eaten as a protein supplement during the summer. 
Irawan et al., (2006) estimated protein content in at 2.23% for D. esculentum. 
However, the present value of (3.84%) was lower than the findings of most of 
the earlier workers viz. Seal (2012), Tag et al., (2014), and Tongco et al., (2014). 

From the above mentioned observations it could be inferred that among the 
different species analysed, D. cochleata was found to have the highest amount 
of crude fat (0.40%) and second highest amount of vitamin C. (23.13 mg g 1%) 
and protein (11.84%). In addition, D. maximum, which had the highest 
amount, could be used as a supplement for vitamin C and T. coadunata for 
protein supplement as it had the highest protein content. Diplazium 
esculentum and D. dilatatum could be used as a fiber supplement. All the 
species had a high amount of moisture content. A dietary food having either a 
high amount of protein or fiber should also have high moisture content. Hence, 
these edible ferns can offer a diet without any digestive problems. 

Dietary minerals.—Analysis of dietary elements revealed that these six 
analyzed ferns can supplement the daily dietary needs of Cu, Fe, I, Mg, Mn, 
Mo, and Zn for human beings. However, these ferns had lower amounts of Ca 
(0.66-1.09 mg g ‘°°) and Na (0.46-3.56 mg g 1%). 

Calcium content.—In the present study the calcium content (0.66—1.09 mg 
g '°°) was low in all the ferns under study. These values were lower than the 
calcium content of all the leafy vegetables under comparison (Table 5). Two 
reasons could be attributed to the low amount of Ca in the analyzed ferns. One, 
the average annual rainfall of Sikkim Himalaya is 3500 mm. Because of this 
heavy rainfall, basic elements, especially Ca, is leached and soils in which the 
ferns are growing are mostly acidic, with low Ca and enhanced Fe and Al 
(Singh, 2013). Also, high Al in soil interferes with the uptake of Ca (Schmitt et 
al., 2017). 

Copper content.—The Cu content in the analyzed ferns were very high, 
ranging from (11.29 to 19.25 mg g °°), whereas most of the leafy vegetables 
ranged from (0.03 to 0.27 mg g '’° ). The unusually high amount was because 
ferns tend to accumulate more Al, As, Cd, Cu, Pb, and Zn (Francis and 
Peterson, 1983; Nishizono et al., 1987; Roccotiello et al, 2015; Van et al., 2006) 
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Iodine content.—Natural vegetable sources of iodine are rare. However, six 
ferns under study were good suppliers of iodine. The amount of iodine varied 
from 0.17 to 0.23 mg g |”, 

Iron content.—The iron content of the ferns showed relatively higher values 
than most of the leafy vegetables such as amaranth (2.32 mg g '°°), asparagus 
(2.14 mg g °°), and drumstick (4.00 mg g '°°). Even other edible ferns such as 
ostrich fern and tree fern had lower amounts of Fe. Diplazium maximum 
(16.35 mg g '°°) was estimated to have maximum iron content followed by 
Dryopteris cochleata (15.96 mg g '°) and Tectaria coadunata (11.01 mg g °°”). 
The high amount of Fe may be because the acidic soils of Sikkim had high 
amounts of Fe due to leaching of bases. High soil Fe could have contributed to 
high leaf Fe. 

Magnesium content.—Among the different fern species, magnesium content 
ranged between 6.54—-13.39 mg g ‘°°. The values could be comparable with 
celery (11 mg g '°°), asparagus (14 mg g '%"), lettuce (13 mg g '°°), and taro 
leaves (5 mg g '°"). 

Manganese content.—The ferns had high manganese content ranging from 
6.54—20.35 mg g ‘°°. Comparatively, the leafy vegetables which ranged from 
0.1mgg ‘°° incelery to 1.06 mg g ‘°” in drumstick leaves, had lower values to 
these ferns. 

Molybdenum content.—Himalayan ferns also showed they are good 
suppliers of Mo, especially D. dilatatum. Compared to common leafy 
vegetables like amaranth (0.13 mg g '°°), coriander leaves (1.12 mg g '%°), 
and spinach (0.01 mg g ‘’’) (Gopalan, Sastri, and Balasubramanian, 1989), 
these ferns supply very high amounts. | 

Zinc content.—The zinc content estimated in all the ferns (0.09-0.23 mg 
~*°°) except D. maximum were similar to most of the leafy vegetables such as 
celery (0.13 mg g‘°’) and lettuce (0.20 mg g '°°). Whereas Diplazium 
maximum contained 2.84 mg g ‘°° of Zn. The unusually high amount can be 
attributed to the accumulation pattern of ferns (Francis and Peterson, 1983; 
Nishizono et al., 1987; Roccotiello et al, 2015; Van et al., 2006). 

Sodium content.—Ferns had very low amounts of Na compared to most of 
the leafy vegetables. Only D. maximum (3.56 mg g '°°) and tree fern had a little 
higher range. This was likely due to the soils in which these Himalayan ferns 
are grown, which are typically acidic and normally had small amounts of Na. 

Heavy metals.—The heavy metals accumulated were compared with the 
CODEX Standard-193 (1995) and it was found that all the ferns accumulated 
higher amounts of As, Cd, and Al than the permissible limit. With the 
exception of D. maximum, the other fern species accumulated lower amounts 
of Hg than the permissible limit. In the case of Pb, three species, viz. D. 
dilatatum, D. maximum, and Dryopteris cochleata were within the permissible 
limit. As the Sikkim Himalyan soils were acidic and rich in Al (Singh, 2013), 
ferns accumulated high Al (Schmitt et al., 2017). In addition, ferns and 
lycophytes tend to accumulate more amount of Al, As, Cd, Cu, Pb, and Zn 
(Francis and Peterson, 1983; Nishizono et al., 1987; Roccotiello et al, 2015; Van 
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et al., 2006). Hence, more accumulation of these elements in a majority of the 
species was expected. 


CONCLUSIONS 


These six edible fern analyzed had better nutritional content than some of 
the commercial vegetables commonly utilized. These edible ferns have great 
potential to overcome nutritional deficiency especially in hilly areas. The 
analysis of dietary nutrients showed that the edible ferns were rich sources of 
protein, vitamin C, and crude fiber. The mineral analysis of these ferns showed 
that they are good suppliers of Cu, Fe, I, Mg, Mn, Mo, and Zn. In addition, if 
their accumulation of As, Al, and Cd could be reduced within the permissible 
limit through cultivation practices, they are safe for consumption. Apart from 
these advantages, they are inexpensive and one or other species could be made 
available throughout the year, if their harvesting season could be extended to 
January and February by cultural practices. Based on the dietary nutrients and 
mineral content, these edible ferns can be popularized and further research 
work could be conducted for standardizing practices so as to exploit them at a 
large commercial scale. In the changing climatic scenario, edible ferns could 
add to the food basket of human beings. 
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Coryphopteris simulata (Thelypteridaceae), A New 
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ABsTRACT.—Previous hypotheses on the evolutionary affinities of the eastern North American fern 
species most recently recognized as Parathelypteris simulata have included hybrid origin from 
North American progenitors and disjunctions with east Asian species. Recent phylogenetic 
evidence indicates that Parathelypteris sensu Ching is not monophyletic and that its species 
belong to two distinct clades, each of which includes members of other widely recognized genera: 
one of these clades includes the type of Parathelypteris, and the second includes southeast Asian 
and Pacific species of the more recently recognized genus Coryphopteris along with some 
temperate Asian species of Parathelypteris. Newly generated rbcL data representing the North 
American taxon Parathelypteris simulata places it in the Coryphopteris clade, as well. In 
consideration of ecological and phylogenetic evidence, the new combination Coryphopteris 
simulata is proposed. 


Key Worps.—Parathelypteris, Thelypteris, disjunction, rhizome morphology 


Relationships among the temperate North American species of Thelypter- 
idaceae and their Asian relatives have puzzled pteridologists for more than a 
century. The Massachusetts fern, Parathelypteris simulata (Davenp.) Holttum, 
is emblematic of such uncertainty, and disparate ideas about its affinities have 
been proposed based on morphological and cytological evidence. 

In the protolog for Parathelypteris simulata (originally described in 
Aspidium), Davenport (1894) considered it morphologically intermediate 
between P. noveboracensis (L.) Ching and Thelypteris palustris Schott, and 
entertained the possibility that it was a hybrid before describing it as a distinct 
species. Wagner (1963) proposed that P. simulata (n=64) indeed was of hybrid 
origin from P. noveboracensis (n=27) X T. palustris (n=35), resulting in an 
allopolyploid with n=62, followed by dysploidization. The hypothesis of an 
allopolyploid hybrid origin between extant North American progenitors, 
however, has not been supported by subsequent data. Smith (1971) instead 
proposed a polyploid origin of P. simulata from an ancestor with base number 
x=32, as then known from the Asian species Coryphopteris pectiniformis (C. 
Chr.) Holttum (Manton, 1955). 

Other ideas about relationships among north temperate Parathelypteris 
also have invoked disjunctions between North America and northeastern 
Asia (see Xiang et al., 2015). These groupings generally ally North American 
taxa with Asian taxa having gradually reduced proximal pinnae, e.g., P. 
nevadensis (Baker) Holttum and P. noveboracensis with P. nipponica 
(Franch. & Sav.) Ching; or ally those that lack gradually reduced proximal 
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pinnae, e.g., P. simulata with P. glanduligera (Kunze) Ching and P. japonica 
(Baker) Ching (e.g., Ching, 1963; Tryon and Tryon, 1974). 

Subsequent phylogenetic analyses (Almeida et al., 2016.; Ebihara, 2011; He 
and Zhang, 2012, N. Patel, S. Fawcett, M. Sundue, and J. Budke, in prep.) have 
recovered two distinct clades containing taxa of Parathelypteris sensu Ching, 
with topologies that do not correspond closely with previously proposed 
disjunctions (Fig. 1). One of these two clades included the type of the genus, P. 
glanduligera, with members of the mostly neotropical genus Amauropelta 
nested within it. The second clade (hereafter, the Coryphopteris clade) 
included all sampled species of Coryphopteris sensu Holttum (1971) along 
with several representatives of Parathelypteris sensu Ching, which in turn was 
sister to the first clade plus members of Metathelypteris, an east Asian genus 
comprising 12 species (PPG I, 2016). 

Based on the most recent molecular phylogenetic sampling (N. Patel, S. 
Fawcett, M. Sundue, and J. Budke, in prep.), Parathelypteris simulata is 
sister to P. nipponica within the Coryphopteris clade, and may represent a 
recent disjunction (Fig. 1). This result is based on the inclusion of two 
widely separated populations of P. simulata, which had identical rbcL 
sequences. Parathelypteris simulata was resolved as the only North 
American member of the Coryphopteris clade; the other North American 
Parathelypteris species, P. nevadensis and P. noveboracensis, were recov- 
ered in the Parathelypteris s.s./Amauropelta clade, in accordance with 
Almeida et al. (2016). This result suggests that lamina shape is not a reliable 
character for distinguishing between these groups, because P. nipponica has 
gradually reduced pinnae similar to members of the Parathelypteris s.s./ 
Amauropelta clade, and in contrast with other members of the Coryphopteris 
clade. 

Throughout its range in North America, Parathelypteris simulata is a plant 
of low nutrient, acid soils, a habitat preference shared with Old World 
Coryphopteris (Holttum, 1982). A notable difference is that P. simulata has a 
slender creeping rhizome, in contrast to the erect caudices in Holttum’s (1976) 
concept of Coryphopteris. Holttum (1976) noted, however, that some of 
Ching’s northern species of Parathelypteris — some of which have subse- 
quently been transferred to Coryphopteris (Zhang, 2012) — had slender 
creeping rhizomes. This pattern of creeping rhizomes in temperate taxa and 
erect caudices in tropical taxa can be observed within other fern genera, for 
example, in Diplaziopsis (including Homalosorus) (Price, 1990; Kato and 
Darnaedi, 1988), and Diplazium (eFloras 2018). Furthermore, Holttum (1976) 
noted that these same northern species of Parathelypteris sensu Ching had 
spore and sporangia morphology like that of Coryphopteris, but recognized 
that “to distinguish them clearly would need more recondite techniques than 
are available to the herbarium botanist” (p. 20). 

Although Holttum made great strides in segregating the tropical members of 
the Coryphopteris clade from Parathélypteris sensu Ching, my ongoing 
phylogenomic research in collaboration with A.R. Smith and E. Sessa, et al., 
seeks to further assess relationships among these taxa. Additional combina- 
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Fic. 1. Bayesian Inference phylogeny based on eleven chloroplast and nuclear loci representing 
203 species of Thelypteridaceae (N. Patel, S. Fawcett, M. Sundue, and J. Budke, in prep.). The tree 
was pruned to demonstrate relationships among members of the Amauropelta, Parathelypteris, 
Metathelypteris and Coryphopteris clades. Posterior probabilities > 90% are shown at 
corresponding nodes. Note that the position of Coryphopteris simulata (in box) within the 
Coryphopteris clade. 


tions in Coryphopteris, still needed for some Asian taxa, await results of these 
efforts. Recent evidence firmly places P. simulata within the Coryphopteris 
clade, however, and forthcoming publication of regional floras in North 
America warrants taxonomic updates. In consideration of ecological and 
phylogenetic evidence, and in the interest of broadening the current 
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circumscription of Coryphopteris to include its sole North American member, I 
propose the following combination: 


Coryphopteris simulata (Davenp.) S.E. Fawc., comb. nov. Aspidium simu- 
latum Davenp., Botanical Gazette 19(12): 495. 1894. Filix simulata 
(Davenp.) Farw., Report of the Michigan Academy of Science, Arts and 
Letters 6: 200. 1904. Thelypteris simulata (Davenp.) Niewl., American 
Midland Naturalist 1: 226. 1910. Parathelypteris simulata (Davenp.) 
Holttum, Webbia 30(1): 195. 1976. Wagneriopteris simulata (Davenp.) A. 
Live & D. Love, Taxon 26: 325. 1977. 
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SHORTER NOTES 


Davip S. Conant, 1949-2018.—David Stoughton Conant was born on June 17, 
1949 to David Perry Conant and Aune Lepisto Conant. He grew up in 
Charlestown, New Hampshire, on the banks of the Connecticut River. The 
Connecticut River runs through David’s life as it runs through the New England 
landscape. As a boy, he spent summers jumping into its swift waters from 
Charlestown’s high bridges. As a young father, he ferried his daughters down its 
winding reaches on a homebuilt raft. And as a professor at Lyndon State College, 
he led students on canoe expeditions to explore its many islands. 

Dave Conant attended Fall Mountain Regional High School followed by the 
University of New Hampshire, where he received his Bachelor of Science 
degree in botany in 1971. At the University of New Hampshire, Dave came 
under the influence of Albion Hodgdon, a plant systematist with a strong 
interest in the alpine flora of New England. Hodgdon was Dave’s style of 
mentor, always ready for a high-risk venture after rare plants and at the same 
time a stickler for high-quality work. It was Hodgdon who encouraged Dave to 
pursue funding for his first field work in Puerto Rico. 

It was also Hodgdon who introduced Dave to his friend Rolla Tryon, a 
prominent scholar of fern systematics and biogeography at Harvard University. 
In 1972, Dave began his Ph.D. research at Harvard under the direction of Rolla 
and his wife Alice, focusing on the ecology and evolution of tree ferns in the 
Caribbean. On completing his Ph.D. in 1975, Dave sought and won a job in the 
Upper Connecticut River Valley at Lyndon State College. With him came his 
partner and soul-mate, Aminta “Minty” Kitfield Conant, whom he met while 
serving as a laboratory teaching assistant at the University of New Hampshire. 
Dave spent his academic career at Lyndon State College, and during this time 
he received multiple National Science Foundation grants. He undertook many 
high-risk ventures to remote and untraveled places to collect specimens, 
including Venezuela, Borneo, New Guinea, Australia, the Dominican Repub- 
lic, Puerto Rico, Costa Rica, Brazil, and Mexico. 

At both Harvard and at Lyndon State, Dave’s intellectual passion was 
developing insights into the ecology and evolution of the scaly tree ferns 
(Cyatheaceae). Dave was a pioneer in applying molecular techniques to 
uncovering the unique evolutionary histories of these tree ferns, focused on the 
plants that he encountered in the course of his extensive field work in Puerto 
Rico. Dave’s approach was broad-based: he made significant research 
contributions in an array of botanical subdisciplines, including reproductive 
biology, ecology, systematics, and phylogeny. Early in his career, he carried 
out what is still one of the few documented surveys of the pattern of spore 
dispersal away from ferns, using radioactive labeling of a Puerto Rican tree fern 
(Conant, 1978). Having documented hybrids in Antillean Cyatheaceae 
(Conant, 1975), Dave went on to propose a completely new mode of speciation 
for the scaly tree ferns, autogamous allohomoploidy, under which new species 
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Photo from David Ballou, Lyndon State College archives 


are produced following hybridization and isolated by self-fertilization (Conant 
and Cooper-Driver, 1980). At the dawn of the era in which molecular data were 
brought to bear on the study of plant systematics and evolution, Dave explored 
the utility of spore storage proteins for interpreting evolutionary divergence in 
tree ferns (Conant and DeMaggio, 1989). At the same time, he substantiated his 
allohomoploid hypothesis with gametophyte culture experiments, demon- 
strating that tree fern hybrids are likely to be selfers (Conant, 1990). The 1990s 
brought chloroplast-DNA restriction-site analysis to the study of fern 
evolution; once again Dave was among the first to explore the new technology, 
applying it to his chosen group, the Cyatheaceae (Stein et. al., 1992; Burke et 
al., 1993; Conant et al., 1994, 1995; Stein et al., 1997). This work provided the 
first molecular data supporting Rolla Tryon’s interpretation, based on petiole- 
scale characters, of some of the early events in the evolution of the scaly tree 
ferns. Major NSF funding made possible Dave’s trip to collect scaly tree ferns 
across several regions of the Old-World tropics, resulting first in a publication 
on the phylogenetic relationships of the tree ferns of Mt. Kinabalu in Borneo 
(Conant, 2001). Then, with the advent of DNA sequencing technology, Dave 
used these new tools to develop a worldwide scaly tree-fern phylogeny, 
teaming up with a broad group of fern specialists to accomplish his goals 
(Korall et al. 2006, 2007; Stein et al., 2010). In one of Dave’s most recent 
publications, he brought his work back home to the Connecticut River Valley, 
contributing a chapter on flora to Rebecca Brown’s Atlas of the Upper 
Connecticut River Watershed (Conant, 2009). 
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As a research mentor, Dave was skillful and generous, welcoming dozens of 
Lyndon students into his laboratory over a quarter century. Many who had 
never imagined themselves as researchers experienced their first success under 
his direction. Some took part in their first scientific meeting traveling with him 
to the New England Fern Conference at Harvard Forest. Dave was also an 
inspirational teacher: his enthusiasm for natural science was infectious and his 
expertise broad. He taught a remarkable array of courses at Lyndon, including 
General Biology, Genetics, The Plant Kingdom, Field Botany and Dendrology, 
Plant Growth and Function, and Advanced Systematics. In later years, he 
continued to develop new courses, such as Deep Impact, a course that focused 
on the catastrophic events that brought to an end the age of the great reptiles. 

Dave’s students looked up to him because he was academically demanding 
and at the same time fair and approachable. Using Vermont’s fields and 
woodlands as his classroom, Dave opened his students’ eyes to the diversity of 
the New England flora. He enticed them with trips to the alpine summit of Mt. 
Washington, the grassy banks of the Connecticut River, and the dripping cliffs 
at Job’s Pond. He instilled in them a sense of wonder and a love for the 
Northeast Kingdom of Vermont: he was a master of place-based education 
before it had a name. Not all Dave’s teaching trips were local, however: some 
were to the fecund wetlands of the Okefenokee Swamp, the shimmering mouth 
of the Suwannee River, and the forested slopes of Puerto Rican mountains. For 
a number of Lyndon State College students, Dave’s field trips gave them their 
first look at the world beyond New England. Those trips transformed their 
perception of the natural world. Many of them went on to become 
environmental, community, and educational leaders in Vermont and nation- 
wide. A number of Dave’s protégées have followed him into the botanical 
profession, bringing his passion for ferns to a new generation of students. 

Dave was also generous in his service to the profession of botany, serving as 
President ofthe New England Botanical Club from 1996—98 and the American Fern 
Society from 2004—06. Both of those organizations prospered under his leadership. 

After retiring from academia, Dave started an excavating business called 
Homework Construction. Dave had a lifelong enthusiasm for heavy equipment 
and enjoyed developing this new skillset in later life. 

Greater even than his devotion to botany has been Dave’s love of the New 
England landscape, and particularly, of his Vermont hill farm located three 
miles from the banks of the Connecticut River. Together Dave and Minty 
stewarded their forests and fields, maintaining a beautiful network of trails and 
wood roads throughout the property. At work on one of those roads, Dave was 
killed in a logging accident on June 27, 2018. He is greatly missed. 
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COVER CAPTION: Asplenium serratifolium Li Bing Zhang & K.W. Xu, a new species de- 
scribed from Bac Huong Hoa Nature Reserve, Quang Tri, central Vietnam. Photo credit: 
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